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ABSTRACT 


The digital computer simulation of the six degrees of 
freedom equations of motion for the XR-3 captive air bubble 
testcraft is presented. The origin of this computer program 
is the SES Loads and Motion Program developed by Oceanics Inc. 
for the Bell 100 ton (100B) surface effect ship. Modifi- 
Cations and procedures used in the revision of the subrou- 
tines to convert the L & M Program from the 100B model to 
the XR-3 model is documented. Measurement data from XR-3 
test runs are used to verify the model for steady-state 
operating conditions in calm water. Computer output for 
turn Maneuvers in calm waters and for regular sea conditions 


are included. 
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ee eRODUC TION 


A. BACKGROUND 

Conventional monohull displacement type ships are drag 
force limited to speeds of approximately 50 knots. A very 
cursory study of the frictional drag component and its di- 
rect dependence on the square of velocity will bear this out. 
It can be shown that the power plant requirements to over- 
come the drag force at the higher speeds requires such large 
plants that the ship weight increases significantly resulting 
in a speed barrier due to drag forces. The desire for speeds 
above the drag barrier speed has led to some new approaches 
Pueship design. The Surface Effects Ship is one of these 
approaches. 

WiEeace Hrrecto. Ship 2S 7a broad Category which includes 
the ACV, Air Cushion Vehicle, GEM, Ground Effects Machine 
and the CAB, oreo’ Air Bubble. Navy Magazine of May 1967 
has a very informative article on different ACV and GEM 
cepts. However, we shall confine our discussion to the 
CAB type of craft. 

Captured Air Bubble Concept - The captured air bubble 
craft consists of rigid side walls with flexible soft seals 
between the water surface and the ship at bow and stern. 

This provides a chamber into which air is blown and cap- 
tured. The resulting pressure in the plenum provides a 
vertical force sufficient to allow shallow draft operation. 

The shallow draft produces a lower drag and a higher 


speed for the same thrust than would be had for a corresponding 
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displacement hull. The reason for this effect is that the 
CAB ship drag is composed of two distinctly different type 
Memeorces The GQrag resulting from skin friction on the 
immersed surface is directly proportional to velocity squared. 
The drag resulting from the friction of the bubble at the 
air-water interface is inversely proportional to velocity to 
the 1.566 power. Thus the sum of these two drag forces has 

a point at which an increase in speed actually results ina 
meauction in drag. 

With the development of faster submarines and the spec- 
tacular increases in speed of aircraft the Navy is rightfully 
Searching for a method of increasing the surface ship speed. 
Ppmeeeceial project office, Surface Effects Ship Project Office 
(SESPO), has been created to specifically look into this 
problem and to evaluate proposals. 

At the present time two 100-ton Captured Air Bubble 
vessels have been constructed. They are the Aerojet General 
100-A and’ the Bell Aerospace Systems 100-B. These two craft 
eC DLOLOLYpeES of a projected POCOstoneeratt.: Pri@r to con; 
struction of the 100-ton crafts, experimental models were 
constructed im the 3-ton range. One such craft, the XR-3 
built by NSRDC, is presently being used in experiments at 
the Naval Postgraduate School. Additionally, Oceanics In- 
corporated was commissioned to develop a digital computer 

\ 
Simulation of the CAB ship loads and motions. The program 


delivered to NPGS in October 1972, was a model of the 100-B. 
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Pee ObJECTIVES 

Eevelopment Oo. the equations Of motion for the CAB 
ship were made in terms of the various elements contributing 
separate forces to the different degrees of freedom. The 
objectives of the simulation given in reference 9 are here 
restated. 

1. To provide a computerized form of six degree of 
freedom equations for SES craft that will yield time domain 
outputs of motions for the on-bubble mode of operation. 

2. To determine the motions of specified, SES confi- 
gurations with special references to human habitability, 
platform requirements (for military or commercial appli- 
cations) and maneuvering requirements for the on-bubble 
mode of operation. 

3. To establish 5 mathematical model for predicting 
maximum wave-induced bending moments, both longitudinal and 
transverse, on SES craft to provide techniques that would 
yield slamming responses for on-bubble mode of operation. 

ve) detects une stability derivatives in the ver- 
tical plane from measured data of the vertical plane motions 
of existing SES craft in a random seaway. 

The objective of this thesis is to closely examine and 
attempt to improve the mathematical model. Our method is 
one of revising and changing the 100-B model L&M program so 
that it accurately describes the motion of the XR-3 model 
test craft. Program conversion consists of subroutine modi- 


fications. In certain subroutines it was necessary only to, 
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input to the program new data, based upon meaSured data from 
the test craft. In other subroutines the method of calcu- 
hatcion, ageb ticle, or entire equations had to be revised. 
The thesis discusses thoroughly these changes and their cor- 
responding rationale. In addition, a comparison of measured 
test craft data with the L&M program output for steady-state, 
calm water operating condition is used to establish the 
accuracy of the simulation. 

A users manual for this Fortran IV program is in Ap- 
pendix A which includes descriptions of the subroutines and 
the method of inputting data. The program is designed to be 
operated in six degrees of freedom. Since many of the 
crafts parameters are input through the use of data cards 
a thorough reading of the subroutine descriptions and their 
ee. eee wiiputs should be made before attempting to 


exercise this program as it presently exists. 





eee eOouUR TIONS OF MOTION 


fen PRODUCTION 

Motion of any vehicle which is free to move in a three 
dimensional space is made up of six degrees of freedom. 

That is, the vehicle (in this case the XR-3) may have trans- 
lational motion in any one or more of three coordinate di- 
rections and it may also have rotational motion about any 

Of these coordinate directions. In order to study these 
motions and the control problems of the XR-3, a mathematical 
model consisting of these six equations is required. These 
equations can then be used as the basis of a computer model 
Suitable for simulation studies. In general, these above 
equations will be complex and nonlinear. The initial de- 
velopment and computerization of the six degree of freedom 
motion equations was based on the use of model test data, 
limited XR-3 fest data, and theory for estimating the co- 
efficients or constituent terms that appear in the equations. 
The mane rent motion characteristics of the XR-3 are determined 
by exercizing this computer model. These results are com- 
pared with available test data and finally used to predict 
other motion characteristics not readily obtainable through 
Mus eCrart Lesting. 

Early analyses of motions within the vertical and la- 
teral planes were considered separately and revealed that 
linear equations provided an acceptable description of 
MOtloncewithin the vertical plane (the symmetrical plane for 


the XR-3 as shown in Figure 1.). However, the linear lateral 
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Figure 1 SES Plane of Symmetry 
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plane equations proved to be inadequate and it was clear 
that nonlinear effects were significant in the description 
of asymmetric motions. 

wOouPLroduce a resultant functional form for the equations 
that will be the most simple form possible in order to achieve 
an ease of computer representation and solution while still 
including the main factors that affect craft motion, a sim- 
plified nonlinear representation was developed. 
COORDINATE SYSTEM 

The equations of motion describe the craft's movement 
in terms of an inertial reference frame in order to determine 
its position in space. Since the XR-3 is expected to have 
small angles of pitch and roll, the main angle indicating 
Orientation of the ship relative to the inertial frame is the 
yaw angle. The range aii this angle is large in order to pro- 
Pemin cgescribe craft turning motion in the horizontal plane. 

A second set of coordinate axes, used to describe craft 
roll, pitch and heave, is the local level reference frame. 
The horizontal ane of this reference frame remains parallel 
to the horizontal inertial reference plane however its origin 
1s placed at ene craft's center of gravity above the mean 
free surface and its x-axis is in the craft's heading di- 
rection. The transformation between the local level and the 


inertial reference planes is given by 


uw ‘ 

te cos¥, Sine o a 
: Ws 

A N 

ky, 0 f 0 f 1 Ky 
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The craft's heave, pitch, and roll motions are related 
to the local level reference plane through a third set of 
coordinate axes called the Body Frame. The Body Frame, also 
called the ship axis system, has its origin at the ship's 
center of gravity and its positive x-axis is also aligned 
in the direction of the ship's bow, the y-axis is positive 
to starboard, and the z-axis is positive downward. The 
transformation between the body axes and the local level 


axes iS given by 


A “ws 
1B Axx Axy Axz ty, | 
A > van 
JIp| = |Ayx Ayy Ayz Wie 
I bn 
Rp ZX Aoy Anz a7 
where: A = cosé 
one 
Ayy = Oo 
ooae = -sin8g 
Rox = sinfsin8 
A = Os 
YY AOE 
A = sin@gcos® 
Wie 
Ax = cosgsin8 
Aly = -sing 
Agog = cos¥%cos®@ 


Finally, to describe the orientation of the fixed body 
frame with respect to the inertial frame a set of Single 
axis rotations are performed. These rotations, through what 


are known as Euler Angles, bring the inertial frame into 
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alignment with the body frame. Taking the sequence to be 
positive rotations about the , 3 and 7 axeSe ln chat sOorder , 
the corresponding rotations areW (yaw), ® (pitch), and §@ (roll), 
and the transformation between reference frames is given in 


term's: of these Euler angles by: 


& an 
1B Byx Byy Byz aT 
- yn 
ae P yx yy yz IT 
k f 
B Box Boy Boz ai 
where: 
eee = cosOcos ¥ 
= cos@sinY? 
xY 
SS ae Sing 
"se = sinfsin@Ocos ¥ - cosfsinY 
Bea sinfsin@sinY + cosfcos V 
Buz = singcos® 
oe = cosf#sinécos ¥+ singsin 
re cosfsin@sin¥Y - sinficosY 
Boz = cosfcosé 


Figure 2 illustrates the relative orientation of the above 
mentioned reference frames. 

During normal operation of an SES vehicle, roll angle @ 
and pitch angle 98 remain within a range where no great error 
is introduced by using linear approximations for the tri- 
gonometric functions. The resultant transformation between 


the inertial frame and the body frame now becomes: 
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Se {ie — 
7 7 a Al 7 . : i 
wey G A ¥=0 LOCAL (TRANSLATED x, 
| | a Lvs v FRAME YAWWED 7) 
Lv= 0 | 
L | 
! 
! 
[ ee ey 
—— 
a 
| } i, | a 
va | es Os Bo 6 
Ue Jp AY 5 ip | BODY (PITCHED @ 
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i. pad . b9= 0 
on ee 
| Peon 
ee 
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a cos, sinY, -0 et nef 
Se = -sinY, cos, g - - 
rS , 

K 8cos¥ + Bsin™, O@sin¥ - cos, el k 


DYNAMICS AND KINEMATICS 

There are three general categories of forces acting on the 
Memo eype Craft. The first category of forces are those acting 
upon the hard structural elements of the ship (rigid body 
forces) such as the sidewalls, rudder, and deck. The second 
group of forces are those acting on the seals and transmitted 
to the ship through their attachments. Finally, the last 
Category are those forces which act on the ship and water 
Surface due to the air bubble and the lift fan system. Figure 
3 iS a Simplified illustration of three force categories acting 
mae obo Lype craft. 

PeuatiemMoOcion OCCUrS aS a result of a combination of all 
three force categories. When the vector sum of all these forces 
are in perience: the craft will continue in a state of uniform 
motion. Motions in the plane of symmetry are: forward and aft 
MOELON OF SUrge; vertical motion or heave; and pitch motion or 
rotation in the vertical plane. Motions in the lateral plane 
are; sway motion or lateral translation; yaw motion a rotation 
Boome tne vertical axis; and roll motion or rotation about the 
forward and aft axis. 

The craft will accelerate in the forward direction as a 
result of propulsive thrust at some equilibrium velocity. The 
Pitch motions of the craft are determined by the interactions 


Sormectmenoust, drag, lift and gravitational forces. Heave 
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motions result from variations in bubble pressure, buoyant 
forces on seals, sidewalls, and waves reacting against the 
weight of the craft. 

Motions in the lateral plane result from sideforces 
acting through the center of gravity to produce lateral 
translation or sway, Sideforce offset longetudinally from 
Pmemcenter Of gravity to produce angular yawing motion, and 
lift force offset laterally as well as sideforce offset 
wereically from the center of gravity to produce angular 
motling motion. 

In order to describe the above motions in a simple and 
accurate manner which is conducive to computer programming, 
emeerrect choice of kinematic variables and coordinate sys- 
tems must be made. When the equations of motion have been 
solved for that set of variables, all other variables can 
then be evaluated Neva ieie relationships between kinematic 
variables; and when the orientation of the body frame has 
been established, all other vectors however expressed can be 
resolved into any common frame using the previously given 
transformations Reser coordinate systems. 

The equations of motion have been based on the following 
chosen kinematic variables: the Euler angle, and their time 
derivatives: wi 9,8; V,8,8; WV, 0,0: and the components of 
the velocity of the origin of the body frame, resolved in the 
direction of the local level frame axes; and the time deri- 
vatives of these velocity components, Uy, Vy. Wy, and 


: ‘ : A p ; 
Uy. Vir Wy, where u; is the component in the i, direction, 
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Ves on the Ir, Glrecrion, and w 


: F “N . : 
L 1s in the Ky aarecevon . 


L 


The equations of motion are to be solved for 


Y, 0, g, ane Vie We 


two advantages to this coice of variables and reference 
memes are; first, the simplicity in the integrations re- 
quired to determine the orientation of the body frame and 
tne location of its origin, and second, is that it. facili- 
tates a reduction of degrees of freedom. That is, the e- 
quation for Z can be eaSily reduced to obtain a representation 
of the heave-only mode which is a readily interpretable 
motion. 

One eeeapabie disadvantage to the above choice is 
that the transformations required to resolve components into 
local level frame axes and to express relationships in terms 
of selected kinematic variables are nonlinear. In most 
cases this process has been simplified by using only the 
imal gsm ee terms in the Euler angles. It is to be noted 
Pia in comparison to the bubble lift force, the other 
forces acting on the ere are generally small; hence, first 
order transformations from body frame to local level frame 
axeS actually provide representative second order accuracy. 

As previously stated a net unbalance of forces or mo- 
ments along or about any axis will cause the craft to ac- 
celerate along or about that axis. This acceleration is 
Galled)a Kinetic reaction. This principle known as D'Alemdert's 


principle can be more formally stated that when the kinetic 
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reaction contributions are included, the sum of all forces 


and moments vanish. Stated mathematically we have: 

oe xk oer, 2 
oy al 
at a: 

A 

0 = Z2-4F, k 
ay 

0 = K-Zm, 7 
asa 

0 = mM-&m. Z 
ieee 

o&~ 

0 = N-2M, k 


where X,Y and Z are the kinetic reactions and BS a F. > 


ah f 


and FP. ‘k and M. 7, M. 5, and M. ¥ are the various forces 
and moments due to thrust, drag, gravity, bubble pressure, 
and buoyancy in the X,Y and Z directions of the local level 
menrerence frames. 

imevaitiatic the Kinetle reactions 1t ks assumed that 
the moments of inertia of the craft are constant with re- 
spect to the Weal level coordinate ssystem= alow, Ee direct 
application of Newton's laws of motion for a rigid body re- 
lative BS ss coordinate system leads to the following 


equations of motion or Kinetic reactions: 


m (u-v¥) = xX 

m (7 + u¥) = Y 

mow = 2 

1 8 - 1,,¥- aE: = XK 
tyy 8 ~ Tye (h ° cS Tex a 
Ive y I., 8 = N 


ae | 





Gemputations Using analog computers have revealed that the 
nonlinear inertial coupling terms included in the above e- 
quations contributed only negligible differences to the re- 
sponses of expected forcing functions (see reference 9). 
Hence, the decision was made to simplify the eguations of 


motions to be used by deleting all nonlinear inertial coupling 


terms. The resulting equations are then given by: 
ma = xX 
mv = OS muty 
mw = 2 
Teg! 7 ro ae 
I 6 = M 
a 
Io¥- 108) Tie 
Note: The u WY product is retained Since it is not a small 


quantity . 


The above simplified equations can be represented in 


° 


Meer ix Lorm by: 


m 0 0 0 0 O u X 
joes. 0 0 oO 0 V Y - mu 
0 0 m 0 0 0 W 7, 
0 0 0 I OP bee g K 
De Oe ae f) 
0 ay M 
foe Oo =r) OUT / N 
XZ Zak 


This matrix representation is used to establish a procedure 


by which the digital computer solves the equations Ge MOuLon. 
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Detailed descriptions of the component forces and mo- 
Ments which contribute to the motion of the XR-3 craft are 
provided in the subroutine descriptions that follow this 


Section. 
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III. PROGRAM SUBROUTINE MODIFICATIONS 


A. MAIN 

The only changes to the main portion of the program was 
the introduction of the term VEL and the corresponding for- 
temenange SO that velocity is written in knots instead of 


feet per second. 


B. SUBROUTINE BOW SEAL 

The XR-3 is presently equipped with soft air-spring type 
bow and stern seals. These seals were designed and cons- 
tructed by NSRDC and installed by NPGS. The seals illustrated 
in Figure 4 are double-celled air bags consisting of a 46" X 
120" frame and fabric ends, sides, and bottom. Twelve e- 
qually spaced steel spring stiffeners reinforce the bottom 
of the seals. The divider membrane has holes in it which 
serve aS a damping device. Winches and cables have been 
attached to the rings slhrwn in Figure 4, so that the shape 
of the inflated seal may be changed. 

To inflate the seals ducted fans are used. The char- 
acteristics of these fans are shown in Figure 5 and the 
Guct ing arrangement in Figure 6. As can be seen from 
Figures 5 and 6, two ducted fans supply air to the bow seal. 
liemrare OL flow is 125 cubic feet per second per fan; hence, 
the total flow to the bow seal is 250 cubic feet per second. 
The return ducting, shown in Figure 6, allows the vane axial 
fans to deliver flow and not cause the seal to inflate like 


a balloon. There is sufficient fan flow to produce a pressure 


30 





yp eangty 





NYSLS PUo MOd 
eign! Laos 
NOILOSS Ws 


oe 





G oznSty 


SOILSIUSLOVUVHD MOTA Nwa 


(DaS/gl4) = MOU «Z9VAVA Awlal « “o 


Sy a Osi Sha) O01 OL Os St Qo 


7s 


SYNSS3Yd 3B89nAa wv 






(AS3) 


WS) (oleae i BEC Ta) isc o) 





SAva AMA WOOF 
(25 S@0Qn DIULS 





NOILI0224 SISLOYUNY aV9NI7 NON 


Fu aaads DWYS LY SSNIONZ NWS 2an Nv 


>SNOILLIQNOD OS 


BZ 





—— ee ee 


| OQ Sas Fa 


NOILVSNDISANOD SONILONG 





"we = eee es ao ee Pe ee och ee ; ? ‘ S Se 
: : — rat t ~ aan =~ Gf ey tam we = ee oe : ? - a 
’ a Voge A, 4 : 
‘ ; .7 a s" *s 

) a ead Ae ES — ew ae ee we ee eee ew ot oe ae ae | 


oe ee eg eee ~ 
2 IC 
eee a Go we enon nn iS 


P 








ee 
apes 
ad ee een tas eee ee Saat haste gma, Oem 7 None 
eh -fa as 

eam on = we ies —- —~ ane t i a See 
222 ee ee ieee 
b ae " J 

DoD Co Coo oot aay eee | - — ies 
. SP aay 

: | 7 Ie 

ee Las ) 

— 8 = wf ow ee we oe ew ee ew ap ee ew ow of a 4” 
=e a Te 


J re ly. 
~ a a: 


33 





differential between the seals and the plenum chamber of 
about 1 to 2 psf. depending on the plenum pressure. 

fveeoee~n seal of the 100B is of similar construction 
to the bow seal presently installed in the XR-3 and there- 
fore the 100B Stern Seal Subroutine was used as the basis 
for computer representation. By comparison of Figures 4 and 
7, the XR-3 seal configuration and the 100B stern seal con- 
Figuration respectively, we can note strong similarities. 
Both are of pressurized fabric-type material that is assumed 
tO have no inertia or dynamics and are constructed with a 
vented membrane separating the lobed air bags. The throttle 
shown in Figure 7 is merely a controllable path for the air 
to return to the plenum. This throttle corresponds to the 
return ducting previously mentioned. 

The amoeeee of Odenton of the program involved 
the identification of the elements of the 100B stern seal 
Calculations so that proper subsititutions could be made 
without destroying the continuity as well as the way in which 
the program coupled with itself. 

The first computation considered was the calculation of 
the bow seal gap. GAP is here defined as that space resul- 
ting from lifting the seal beyond the surface of the water, 
Miatets, an alr gap resulting in a loss of air from the 
plenum and increased air flow. The equation is here repro-~- 


duced without its iterative subscripts. 
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GAP = -(Z+ZS-XX (3,K) *THETA+YY (3,K) *PHI+ETA (3,K) ) 
XX (3,K) = XBS+RBS*COS (ANG) 


Mayet 3, K) 


RBS* SIN (ANG) 
fieeocalculation results in the distance DS as shown in 
Figure 9. 

dhe purpose of the XX(3,K) and YY(3,K) calculations are 
to locate points along the perimeter of the hinge. The dis- 
meamecesS XBS, XX(3,K) and YY(3,K) can all be found in Sub- 
routine INCON. RBS is not required for the XR-3 seal and 
is removed. The reason for the calculation of RBS was that 
the 100B hinge is semicircular in shape as shown in Figure 8. 
The XR-3, on the other hand, has a seal that has its hinge 
On the bow of the ship and is affixed straight across, ath- 
Wartships, therefore XX(3,J) = XBS. The distance, YY(3,Jd) 
is an mer enental aaeueriee athwartships and is also cal- 
culated in INCON. 

YY(3,d) = 0.5*XBBW+(J-1) *DELYBS 

In this way we have located a series of straight line seg- 
ments running athwartships at the seal hinge location. 

“eallidediiy the overhang of the bag portion of the 100B 
bowseal had been calculated as ELBAG by Subroutine FGl. 
Since the XR-3 bow seal has no overhang, this calculation 
has been removed. The bow and Sears, seal of the XR-3 are 
basically alike and therefore the only @itterences are in 
Piemoietinition of terms used in the equations. The cor- 
responding changes have been made in INCON, RUDDER, RHS, 


ana wherever else the Common Block ties were made. New 


AG 





ES es eee OE ee ee 


ee —<-_ "Se 2 Ca 6 Ee ewe e+ Oren fe eee ee Se 


[ 
| 
[ 





lOO 6 


Figure 8 





6 9anbty 


NOILVYNSISNOD IV3S MOE 


O'O = W104 : 
O'O = HO1Id MBIA Jails qais 
WYVW LIVYO 99 
8 GNV 2 SNOILVNOF Ag a i 


Q31V901 S! LNIOd SIHL a 





NICaeao NIG iy es 


x8) 


Y3IGNVHO 

WAN3 1d —s 

JO dol 

aM4 ¢——— 
$—___—_—__——- Sax : SX — 
| 
ISEHL 
JIONY vas 


WOSNVYL OL ISEX 











terms that were included match their corresponding terms in 
the 100B Stern Seal Subroutine. With reference to Figure 9, 
theses are XBF which is the length of the seal from hinge 
memeransom, THBSI which is the angle between the craft ver- 
tical and the seal, BLEAK which is the seal base leakage area, 
ZBS which is the distance from the CG to the seal hinge, 

COSBS and SINBS, which correspond to SINTH and COSTH in the 
100B program. These later two terms are the SIN and COS of 
the THBSI bow seal. 

The following is an equation by equation discussion of 
the XR-3 Bow Seal Subroutine. The first few statements serve 
memzero the force, GAP vector, SKI vector and leakage area. 

(1) DELPBG = PBS — PB 

Equation (1) calculates the differential in pressure 
between the bubble Roenin PB, and the bow seal bags, PBS. 
PB is calcualted in Subroutine RHS and PBS 1s an INCON Sub- 


memcane input. 


(2) PBAR EB = sear 


(3) DELP PBAR 


(4) IF (DELP.LT.0.0.) DELP = 0.0 

Equations (2) through (4) calculate the pressure dif- 
ference between the bubble, PB, and the sea level value PINF, 
insuring that DELP is never negative since we cannot create 
a vacuum in the plenum. 

(5) SINDIF = SINBS - COSBS * THETA 

(6) COSDIF = COSBS ~- SINBS * THETA 


COSBS and SINBS are calculated in Subroutine INCON. 
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‘THETA is the pitch angle and is considered small enough to 
use the small angle approximation for the sin. Therefore, 
equations (5) and (6) represent adjustments to the INCON 


meaues, COSBS and SINBS, due to pitch angle. 


{ce 


(7) X1 XBS+2ZBS* THETA-XBF*SINDIF 


(8) Z1 ele Fboraso* Pie TA-KBF*COSDIF 


Equations (7) and (8) heGaneram line on the water ™s sur> 
maee OL the seal based upon the draft and pitch angle. Re- 
ferencing Figure 9 which assumes no pitch or roll it can be 
seen that the distance determined by equation (7) is indeed 
X1. Similarly we see that equation (8) locates the distance 
Zl. It is assumed that the forces acting upon the seal and 
the craft act along this line which is centered in the wetted 
area at the water's surface. The assumption is that the seal, 
being a’ buoyant body Aeee not penetrate sufficiently the 
water's surface, but floats upon it. 

The next task is to calculate the wetted area of the 
seal and the gap, if any. The seal may be divided into as 
many as eleven Peers eecereing to NSTA (3) . 

(9) ELSKI (K) = (ETA(3,K)-DETABX (K) ) * (XX(3,K)-X1)-Z1+ 

VN ( SS Ceo AB IE 

(10) GAP(K) = ELSKI(K) 

(11) IF (GAP(K) .LT.0.0) GAP(K)=0.0 

ELSKI represents the wetted length of the small segment 
SBerrectea for waves, (ETA(3,K)-DETABX(K)): and roll, PHI. 
Equation (11) provides that as long as ELSKI is negative no 
Gap exists. Equation (12) computes an average of ELSKI over 


the segments previously determined. 
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SeaepesokiA = (BEESKE (ita) +ELSKIA(J)) /2.0 

Obviously if the average wetted length computes to be 
less than zero there is no need to make any of the following 
computations and gap is the only meaningful term. However, 
as is usually the case, ELSKIA has a value and it is then 
necessary to next test ELSKIA and determine if it exceeds 
the maximum wetted length permissable. If so, it is then 
Set to the maximum value. This is reasonable since the 
Maximum wetted area would exist when the top of the plenum 
was contacted by the water and the computations would stop. 
The succeeding equations calculate a coefficient of drag 
based upon the craft velocity and then the resulting drag 
forces. Similarly the remaining forces and moments are 
Calculated. 

Mice leakage flow, om was not permitted in the 100B 
fombess the seal broached the surface and GAP became non-zero. 
fmecem@asce or the XR-3 it was found that there existed a 
leakage area between the sidewalls and the bow seal itself. 
This gap was calculated ence nolics as CEES in Block 6 of 
Subroutine INCON. This small area accounts for a flow of 
Meetesimately 10 cubic feet per second. 

Inputs to Subroutine BOWSL will be discussed in the 


section covering Subroutine INCON. 


Seo UBROUTINE COLFIL 
ig@eweoLt?! Subroutine controls output, both form and 


content. Since substantial changes to the Oceanics L&M 


Al 











Peegram Nave been Made, a user's manual for this subroutine 
is included in Appendix A. The modifications herein des- 
Ccribed were made for the purpose of arranging the output 
in a more convenient form for the purposes of this report. 
Presently time histories of 16 variables may be printed in 
two tabled summaries of values and up to 10 graphs plotted 
in each run. Additional variables may be output by a single 
change to the present format. The graphs may be drawn on 
either the off-line (Galeompi= Or on-line (print) plot. 
Changing this routine necessitated minor changes in both the 
feo and INCON Subroutines. 

The changes in RHS consisted of modifying the method 
of writing the variables on the temporary storage device. 
Only one temporary storage device is used. This was done 
in order that the age cau be searched by variable number 
instead of searching by variable name. 

The changes in INCON necessitated the creation of new 
mmeckS 20,21, and 22. Information input and the form of 
that information oe be found in the user's manual, Appendix 


A. 


D. SUBROUTINE INCON 

tic —INEON Subroutine is the most difficult of all the 
subroutines to discuss since it reads and computes the ini- 
tial conditions of all variables and constants which affects 
how the entire program functions. The thorough understanding 


of this subroutines! function and operation are a prerequisite 
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to exercising the program. To assist in this understanding 
a description of the operation of the Subroutine and where 

necessary full description of the important functional com- 
putations is presented in this séction. 

The INCON Subroutine is divided into numbered blocks 
which are controlled by the read statement number 10. (see 
mependix B). Notice that the MeenOdelkS ta read Qa 3 digit 
integer followed by a 2 digit integer. The first 3 digits 
are called the block number and the next 2 digits the option 
teg. By executing the computed go to statement according to 
the block number the statement next executed can be controlled. 
Consider as an example the number 01601. The block number 
is 016 and the option tag is 01. Hence, the block number 
will cause the program to flow to statement number 1600 by 
the GO TO control and nes be routed.. to statement 1605 by 
Seieton tag number in another GO TO control. Using this flow 
Bemerol all of the necessary information is input to the 
Reegram from the card reader and disk. 

What follows ne Nneiwd leewabyeblocksdescription of the 
impute, but rather a description of those places in which 
changes ee SBiseen made and in which particular calculations 
Smmincerest occur. A block by block explanation of the input 
1S contained in the user's manual Appendix A. 

EmocK 2 Option 1 

This input block is the summary mass and inertia input. 

Option 2 allows a discrete mass input to be read. It is 


important to point out the source of the numbers presently 
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being used for the XR-3. As a first approximation the moments 


of inertia found on page 1 of reference 2 were scaled up to 

the present XR-3 configuration. The craft weight given in 

reference 2 was 4192 pounds and a ratio... of 5813 = 1.39 was 
Aa 


used to compute the following inertias. 


Z 
mx = 294) slug ft 


TY = 9320 slug ft° 


IZ 2 


i300 siug £t 
However, the numbers presently employed are those pro- 
vided by NSRDC and are: 


IX = 2870 slug £t7 


fh="9520 slug £t- 


LZ e 


mvS60 Slug £t 
As can be seen, the difference 1S Minimal. No information 
was Mibiable for the eredice inertia I. and therefore it was 
assumed to be zero value. 
feck 4 option 2 

The appendages section of the program has been removed 
Bemce there are no fins or other appendages on the XR-3. 
Block 5 

The measured data inputs to the stern seal block con- 
Sists of the X and Z coordinate of the seal hinge, the seal 
angle, the length of the seal leading edge, the pressure 
differential between the seal bag and the plenum chamber, 
seal leakage orifice coeffficient and the base leakage area. 

It was assumed that the leakage orifice of the XR-3 had 


mvewscame characteristics as the stern seal orifice of the 
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100B. Thus the Seal leakage orifice coefficient remains 
unchanged. 

Tne base leakage area refers to that area over which 
the leakage coefficient exerts influence. This is the gap 
that allows a flow from the plenum to the atmosphere. The 
leakage flow is the control over the plenum presSure, fan 
Eiow rate and thence the fan horsepower:.calculations. The 
leakage area is distributed throughout the craft and it is 
necessary that this number be adjusted to: bring the above 
three calculations into balance at a steady state operating 
condition. If the base leakage area is reduced it causes 
mmemeratt to either rise, or the pressure to increase, or 
the fan flow rate to be reduced. The number presently being 
input is the result of such an adjustment and the three cal- 
culations see an ee varent with experimental data concerning 
draft and plenum pressure. 

PeLOock 6 

These bow seal inputs correspond to those of the stern 

seal ‘La EEOC 5's 


beee@ck 7 option 2 


The critical Froude Number was obtained as follows. 
is! 


is Se ae ms 
Where: V = speed in feet per second (critical or hump 
speed) 
L = length of the plenum (nominal) in feet 
G = gravitation constant 
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Peock 8 option 1 

Those inputs that correspond to thrust have been re- 
moved and are found in block 16. 
eo@ck 9 option 1 

Although there is no rudder on the XR-3, it was con- 
Sidered necessary to provide information for the rudder 
eeeect produced by the propulsion system. Refer to Figure 10 
which shows a crossectional view of one of the outboard motors. 
Dimensions in the figure are inches and must be converted to 
feet and square feet for input. | 

The first assumption made was that the effect of the 
outboard motors could be separated into two parts, first, 
the forces contributed only by rudder action and secondly the 
forces which result from the thrust. 

The second Saunt ian was that the dashed line represents 
the water surface under most operating conditions. That is, 
at the stern, the water level doesn't Change mueh= at lease, 
not enough to reflect great rudder force changes. Addi- 
tionally, waves do not interrupt flow at the rudder. There- 
fore, all that is required is to determine the respective 
areas and lengths, convert them to feet and enter them pro- 
perly through the card reader. The wetted surface rudder 
meet area, therefore, is that region which is cross hatched 
me Figure 10. 

Rudder span is defined as the vertical length of the 
rudder. This distance is the wetted surface distance. Based 


eon the foregoing assumptions only 2/3 of the block marked 
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A in Figure 10 through which the water line passes con- 
tributes to wetted area. 
Therefore, rudder span = 2 (44) +114 =) ose = (1.21 £t.). 
The rudder aspect ratio is defined as the rudder span 
Squared divided by the area, a dimensionless term. 


2 
1 2 ae es’ 
.68 
Total rudder area is the summation of the cross-hatched 


area; A and B. 


The area of B is taken to be the average width 5.75 
inces multiplied by the length. The area of Block A is 
teem tO be a rectangle 11.75 inches by 4.5 inches with a 
small triangle removed. It is felt that these approxi- 


mations will not lead to any gross force inaccuracies. 


R, = 31.94+66.1 = 98.0 inches’ (.68 feet?) 


Block 10 

The aerodynamic inputs are the reference length and 
referenc width. These dimensions are used to compute a 
reference area to be used in conjunction with the parti- 
cular values of the coefficients obtained from the wind 
tunnel tests. 
eereck 15 

The shears and moments calculations have been removed 
and replaced by a dummy subroutine SAM, since no stress 


analysis will be performed on the XR-3. This removal 
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meovedes a Convenience for the purposes of this report. 
Block 16 ppiion LAG 1 

Entries input in block 16 option tag 1 permit the user 
to input a predetermined time history of thrust on the star- 
board screw aS a map. That is,either constant or variable 
thrust may be selected. The user must input data points 
from the map of the variable thruster and the function sub- 
program FGl is used to determine the functional value of 
Eheust at any given time. A maximum of 24 data points may 
be used to describe the map. 
Pmeer 16 option tag 2 

Exactly the same aS option tag 1, except that these 
entries apply to the port screw. 
Bwock 16 option tag 3 

The eee iielOChmlbGwOoDtlOngtag Geallow the user 
to ine a rudder motion Mate Oottive rudder angles are 
right rudder angles which cause the craft to turn to the 
meng ihit . 
Block 17 

Ontion bags and) 2 permit the user to vary bow and 
stern seal differential pressures as the program run pro- 
gresses. The differential pressure is taken to be that 
between the seal and the plenum. ‘This represents a revision 
EOweEne Original program. 
Biteek 19 option tags 1,2, and 3 

These inputs allow the user to input fan map data and 


mime ontrol air flow in both the seals and plenun. 
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fe sUSROUTINE PROP 

The forces and moments associated with the propulsion 
and steering systems for CAB craft are somewhat unique for 
each propeller and rudder system type. In the original 
Oceanics' program, written to simulate the Bell 100 ton 
SES, the system described consisted of two fixed super- 
Cavitating propellers each with a rudder mounted directly 
astern. Both the propellers and rudders were located close 
to the transom of the craft. In the XR-3, the propulsion 
system consists of two forty horsepower outboard motors each 
is fitted with a two bladed screw. The port screw has a 
right hand rotation and the starboard has a left hand rotation. 
With this system steering is accomplished by rotating these 
motors in the horizontal plane.wwin tcuri, the main thrust 
vector produced by each propeller is rotated in the hori- 
zontal plane to provide a desired turning moment. Additional 
turning moment is provided by the motor/propeller housings 
which are essentially Shaped like a rudder foil. However 
the surface areas of these housings are significantly re- 
duced from typical rudder areas associated with a fixed 
Brepeller system. The rudder characteristics are further 
described in the RUDDER Subroutine. 

iiewsOtacion OL the propeller thrust vector in the 
Mer~rzOntal plane adds a complication factor to the XR-3 
description which was not present in the 100B craft simu- 
lation. Hence, the original PROP Subroutine required major 


eiamges in order to make it properly describe the XR-3 craft. 
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Moeene XR-S) Gescription, aS in the 100B, specific pro- 
peller characteristics such as; RPM, Blade pitch, and the 
various propeller coefficients are not included to determine 
the thrust and torque produced by each propeller. In both 
Simulations the thrusts are either input values or cal- 
culated from the equations of motion for the associated 
input velocities. 

The propeller representation allows either propeller 
to provide any desired thrust and side thrust independent 
of the opposite propeller. This feature 1s available 
through the use of function routine FGl. FG1l is used to 
describe the thrust of each “propeller as a function of time. 
When the program operates in this mode the velocities and 
accelerations of the XR-3 are obtained from the equations 
of motion as 5 result of the input propeller thrusts. An- 
piner operating mode provides for velocity to be input and 
the resultant total propeller Eich ls 'Obcained from the 
motion equations. In this mode of operation, it is assumed 
that the total thrust 1s divided equally between both pro- 
pellers, and function FGl is not used since the propeller 
thrusts are not input. 

The above two modes of operation allow for maximum 
Program flexibility and simplicity. Test Gineaditasar bolt la 
thrust and velocity is available from the installed XR-3 
instrumentation. Thus, the simulation can easily duplicate 


any given set of operating conditions. 
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The force and moment equations for the propeller 
system of the XR-3 are derived from the propeller thrust 
Giagram given in Figure 1l. The figure displays the main 
thrust and side thrust vectors for both propellers positioned 
with right rudder angle a, positive pitch angle Q, and ne- 
gative roll angle ® present. The thrust vectors act through 
em@es port and starboard propeller POctE1TONS >) XD, yD; and zp, 
trom the craft center of gravity in the Body Reference co- 
ordinate system. These thrusts are then converted to forces 
in the local level reference coordinate system for the 
@eeeulation Of craft pitch, roll, and heave and forces in 
meee inertial Reference frame for calculation of craft tra- 
jJectory. That is, the craft yaw angleW, and x and y dis- 
placements. 

The forces associated with the propeller arrangement 
ena Figure ll are depicted vectorially in Figure 12 in the 
Local Level Coordinate system. Employing the simpli- 
fications’ used in the derivation of the equations of motion 
namely that the pach angle @ and roll angle @ are limited 
to small VobGesrm Suc that Sine-0, coss=1, and sin@=9, cos=1 


the following force equations result: 


T12 


= THRUST 1 * COSa 
T2X = THRUST 2 * COSa 
STIX = STHRUST 1 * SINa 
ST2X = STHRUST 2 * SINa 
Tly = THRUST 1 * SINa 
T2Y = THRUST 2 * SINa 
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Subil se 


SLANE 


eZ, 


EZ 


ST1Z 


Sle 


= STHRUST 1 * COSa 

= STHRUST 2 * COSa 

= THRUST 1 * COSa * 6 
= THRUST 2 * COSa * 6 
= STHRUST 1 * COSa * 9 


=o RUar 2 * COSa * 9 


By summation of forces in the three coordinate direc- 


memos the following port and starboard propeller and total 


force equations result: 


FXS 


PAP 


fe x 


ENS 


BYP 


FY 
_FZS 
"FZP 


FZ 


eet | SX 


wie ok 


=i homo nal? 


ae aero LY, 


=O Y to ZyY 


= PYSr+ PVP 


Sea Sle, 
=m 2 t+ ST28 


= FZS + FZP 


The moment equations associated with the preceeding force 


equations are 


FKS 
eRe 

bls 
FMS 


FMP 


given as: 

—wijoe Yeer-erYsS * ZP 
= PZPee9_YeP —- FYP *. 2ZP 
= FKS + FKP 

= -FZS * XP + FXS * ZP 


=erZe >) XP + EXP * ZP 
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FM = FMS + FMP 
ENS = —=PxXS * YP + FYS * YP 
PNP axe = YP -payP *~ XP 


bile PENS +) NP 


Throughout the preceeding representation it 1s assumed 
that the magnitudes of the main and side thrust vectors 
remain independent of the various rudder deflection angles a. 
Mees alSO assumed that both port and starboard propellers 


rotate together through the same rudder defelction angles. 


Pee sUBROUTINE RHS 

Subroutine RHS contains the expressions of the right 
hand side of the force and moment differential equations. 
These forces are calculated and summed in this subroutine. 

Changes in this ayeece cles ConstSt Of Modification of 
the Bubble Drag Coefficient of Friction, Common Block ties, 
Pen center of Pressure Location. The rationale for these 
changes iS presented in the section concerning the presen- 


mation Of data. 


G. SUBROUTINE RUDDER 

The rudder subroutine used in the XR-3 simulation re- 
Mains essentially unchanged with respect to its force and 
moment equations. However, some changes are introduced in 
the calculation of rudder deflection angle ‘a' and in the 
input data to reflect the proper XR-3 mudder location and 


aadmensions. 
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Rudder deflection angle can either be input as a con- 
Stant angle in the input data set or it can be made to vary 
with time by using function FGl. Note, this rudder de- 
flection angle is the same angle input to PROP subroutine 
tO control the main thrust vector orientation. 

As previously noted in the description of subroutine 
PROP, the motor/propeller housings constitute the XR-3 
rudders. Physical measurements of these houSings were used 
as rudder dimensions in the input data set. Each rudder/ 
housing is located directly astern of itS associated side- 
Wall. They are vertically mounted and rotate in the hori- 
zontal plane. 

mie analysis Of the AR-3 rudders assumes that they pro- 
mace Only lateral and drag forces. Lift or % direction 
forces are eonmiaeced negligible and hence they are set 
equal to zero. 

The -development of the force equations associated with 
the rudders of the XR-3 follows the analyses from references 
9 and 7. The lateral force produced by both rudders is 
Been by 

Eee AO UU er (1+d/h) *RCLB* (a-(1+d/h) *V/U) 
where A. is the rudder area, RHO is the water mass density, 
Mes the velocity in the X direction, dis the draft at the 
fader, h is the draft plus the rudder span or length, RCLB 
is the lift coefficient for the rudder given as 
RCLB = 2AR/(AR+3) where AR is the aspect ratio, a is the 


rudder deflection angle, V),, 1s the first order approximation 
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mudder velocity in the lateral direction Vy 1 See Wem oy ; 


VT VtXAR*R-ZR*P 


with XR and ZR the coordinate locations of the rudder center 
with respect to the craft center of gravity. 


The drag force produced by the rudders is given by 


FX = -2*CD*AY + RHO*UZ. 


mae drag coefficient is given by 


CD = 2*CFR +PI/8 *(Tx/C)2 * (14G(h-a)/U2) + 
RCLB*V, /U 


where G is the gravity, Tr is the rudder thickness, C is 
mgiemrucder cherd, and CFR is the coefficient of friction 


given by the Schultz-Grtinow formula 


eHRu= 4077 (hog REY — .407)-° °° 


where REY'= U*Ar/(h-d) /ENU. 
REY iS commonly referred to as the Reynold's Number and ENU 
is the kinetic viscosity Ase 10° -£t*/sec. 


The two preceeding force equations produce the following 


MmOMents : 
FK = -Z2R*FY 
FM = ZR*FX 
Pie RAEY 
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H. SUBROUTINE SIDEWALL 

Changes made to the sidewall subroutine consist of the 
removal of force computations associated with ventral fins 
which are not present in the XR-3 testcraft. The remainder 
of the subroutine is unchanged from the Bell 100 ton craft 
version described in reference 9 

The hydrodynamic and hydrostatic forces and moments 
acting on the sidewalls are found by application of slender 
body theory. This theory, based on motion of a foil ina 
viscous medium, developes methods for finding the lift and 
drag forces associated with moving underwater bodies. These 
methods are primarily linear in nature as shown in refer- 
ences 1 and9with the effect of instantaneous changes in im- 
meérsed sectional area, draft, etc. due to orientation changes 
arising from heave, pitch, and roll motions reflecting the 
most significant nonlinear variations. The sidewall forces 
peesormcOontaln certain nonlinear terms due to cross-flow drag, 
which are'‘important for the case of very low aspect ratio 
lifting surfaces ee the XR-3 sidewalls. The sidewall 
forces will then be composed of two separate contributions; 
first the gacne due to sidewall buoyancy and slender body 
hydrodynamic reactions and second the effect of cross-flow 
mag, 1.e. 

Ie = FH see )D, 
Sidewalls SW Sw 

where FH forces result from slender body theory, and FD forces 


are the result of the cross-flow drag terms. 
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Meeieugi me Changes were made to the sidewall sub- 
routine force and moment equations they were considered 
of adequate importance to be hére’ discussed. 

Before computing the forces on the sidewalls, the 
Bpeoroutine first calculates the draft of or the gap be- 
neath the sidewalls. The draft/gap is found for each X- 
@reetion section after it has been corrected for craft 
roll and pitch angles and waves. If gaps are present, then 
the total sidewall leakage area ALSW is computed by summing 
the products of gap heights and sidewall section lengths 
fen both sidewalls. 

1. Cross-Flow Drag Terms 

The cross-flow drag force components are assumed 
Bemoniy be significant for the lateral force, yaw moment, 
and roll moment. While such forces are also present in the 
vertical force and pitch moment, they are not as large re- 
ierve to the predominant forces due to bubble pressure and 
are neglected. The lateral force due to cross-flow drag on 


fecangle sidewall is given by: 
oa ie 2 ~CHOC*YREL*ABo (VREL) *DELX* DOWAV 


where VREL = V + XAVG*R -((ZS-DSWAV) /2) *P and DELX is 
the length of each sidewall section, DSWAV is the draft of 
the ith sidewall section corrected for waves,(ZS-DSWAV) /2 
is the average Z distance of the ith sidewall section from 


XREF where XREF is input as 9 feet forward of the transom. 
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These expressions are appropriate to the particular 
Values associated with sidewall sections in either the port 
Or Starboard side. The yaw moment due to cross-flow drag 
is 


FND = FYD *XAVG 


where XAVG is the average X distance of the ith sidewall 


ewement form XREF. The roll moment is given by 
Pips — FYD “(2S - DSWAV (7/2) 


Tre total force and moment values due to cross-flow drag 


are then 
EXD o= “FYD Clea Fy Oot 2) 
FND = FND ‘(1) + FND (2) 
FKD = FKD (1) + FKD (2) 


where the subscripts (2) and (1) indicate the port and-star- 
board sidewalls respectively. 

The cross-flow drag coefficient is the value cor- 
m=e<ePonaing to that for flow normal to a long flat plate, 
Pith a Vole sor CDC} = 1.28. 

Zeeeolender Body Theory Terms 

Following the procedures indicated in reference 9, 
together with evaluation of all integral terms involving 
derivatives of sectional properties by means CmetMced rat lon 
by parts the following results are found for a single side- 
Gael {port or starboard, using appropriate values for each): 


The X-direction sidewall drag force is given by 


om 





PX kney 2°eDT SawAREA * U * U 


where WAREA is the wetted area of a sidewall section cor- 
Meetea for draft, roll, and pitch. CDT is the computed 
drag coefficient. 

The Z-direction sidewall force corrected for added 


mass effect is given by 


Pane —-G beOe-l*UMASSS* THETA - U *A33S *W + Q * 


UR mec om meso ——8262) —- USAS3S * P *YLSW 


where BCO is the mass of the water displaced by the ith 
Sidewall. element which causes a negative buoyant force, 
A33S is the vertical added mass of the ith sidewall element 
Given by A33S = RHO * TT* BS **2/8 , where BS is its beam. 
BC2 is the total sidewall vertical added mass and YLSW is 
the lateral distance of the ith sidewall] element form the 
Centerline. 


The hydrodynamic lateral force is given by 
FYH = “RODE *U * (Vv + XSS * R- ZS * P ) 
where A22S is the lateral added mass at the stern given by 
A22S = RHO* .4 * PI * DSS**2/2 


migere DSS is the sidewall draft. The moment about the Y-axis 
is given by 
FMH = -U *XSS*XSS* ASS3 * Q + G*BCOO + U * 


(A33 * XSS+BC2) * (W+U *THETA+YLSW*P) 
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where BCOO is associated with the sidewall volume between 
meee and XS. This computation is available by entering a 
table of values computed offline in program SIDETAB and 
stored in common block WAVTAB by the INCON subroutine. 


The moment about the Z axis 1S given by 
FNH = FYH *XSS -U * ((V-ZS *P) * BCS + BC 6*R) 


where BC5 is the total lateral added mass and BC 6 is the 
moment associated with the total lateral added mass. 
The total slender body theory and crossflow forces, 


and moments for both sidewalls are then given by 


FX = FXH (1) + FXH (2) 

FY = FYH (1) + FYH(2) - FYD 

FZ = FZH (1)2+ FZH (2) 

FK = (FZH (2) - FZH (1))* YSW - ZS * FY + FKD 


eh 


(FMH (1) + FMH (2))+4 2S * FX 


PN ee Glee) (2) *Yow + FNH (1) + FNH (2) +FND 


fetne calculation Of the individual sidewall force = 
terms it is understood that the computation iS appropriate 
to the instantaneous values of the separate elements as the 


vehicle moves. 





fr oUSROUTINE STERN SEAL 
Since the bow and stern seals of the XR-3 are similar 
except for location the subroutines also are similar except 
for variable names. 
Pemesence tO Pagure £3 and the list below will acquaint 
the reader with the important variable names. 
Dimension Ohargin 
AS ENCON ROUTINE BLOCK 225. XS is calculated 
by locating the CG and there calculating 
the distance to the transom. 
2.Geveyl INCh ROUTINE "BG0CK 500. This 1S input as 
the distance from the transom to the hinge 
point of the seal. 
ZS LiCC wnOULINE BeOCK 200. 2S is calculated 
in the saine manner as XS except that it 


measures the vertical distance from CG 
eOnecne. Keel. 


4551. .« INCON ROUTINE BLOCK 500. This is the ver- 
Piet Tdictance tromucne hinge to the keel. 


XSS This is calculated in INCON BLOCK 1302, 
using the information XSSI and XS and 
represents the horizontal distance from 
the center of gravity to the seal hinge. 

ZSS Pines tsecalculated in ZINCON BLOCK 1302 using 
Mice IiternatLon provided by ZSSI and ZS 


and represents the distance from the center 
of gravity to the hinge. 


MS oul This is the angle of the seal face made 
With, the Craft vertical and is input in 
block 500 of INCON. 
Lift and drag forces in the stern seal result from the 
pressure differential between the seal and the plenum and 
the amount of the seal in contact with the water respectively. 


Since the plenum exhausts mainly under the stern seal we 


assume the bubble to extend to the tip of the stern seal. 
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Thus the stern seal becomes an air lubricated surface and 
the pienum pressure 1s extended all the way under the seal 
to the atmosphere. 

We can then assume that the forces acting upward con- 
centrate themselves at the point (B) in Figure 13 which re- 
Beets the mid-point of the wetted surface. It is also 
logical to exert the drag forces at this same point, Since 
they act along the wetted surface, the surface of friction. 

We find that the forces in the Z direction are much 
less than in the bow seal because the differential pressure 
is so much less than in the bow seal. The bow seal iS acting 
against atmospheric pressure while the stern seal is acting 


against the plenum pressure. 


See COMMON BLOCK TIE MAP 

The Common Block Tie Map shows which common blocks 
appear in each subroutine. 

The map tS seod ey £inding the common block in the left 
hand column and thén reading to the right. An X appears in 
each column under the subroutine in which that common block 
appears. 

Some common ties have only one entry. This results 


meom the removal of the SAM subroutine. 
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COMMON BLOCK TIE MAP 
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SUBROUTINE 
NAME 
BLOCK NAME 
AIR 
AXIS 
ipMiCO 
COLUMN 
CONST 
eNTRL 
CURVE 
ENGINE 
EQNCO 
pee RO 
FAIR 
FANMAP 
MORBS 
FORSS 
PE ROP 
FROUDE 
E RUD 
GBOW 
GEOM 
GEOMBS 
GEOMSS 
GEOMSW 


COMMON 


oy 





COMMON BLOCK TIE MAP (CONT.) 
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SUBROUTINE 


~*~ ~~ < 
GG tA 
CeCEEe er 
snoop | [| 
ed PLD Pee rer 
COC EC EET 
CCPC EEC 


NAME 
BLOCK NAME 

GRAF 
HEADG 
KSWITCH 
LEAKER 
MASSES 
MATRIX 
MSIDW 
MWAVE 
OPTION 
PLENUM 
PRIME 
PRTINT 
PWAVE _ 
RISER ' 
ROLL 
RUDDR 
SIDE 
SOF TBS 
SOFTSS 
STABLE 
STSLR 
SUM 


COMMON 
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COMMON BLOCK TLE MAP (CONT.) 
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NAME 
COMMON 
BLOCK NAME 
VALOLD 
VARBLE 
WAVE 
WAVEF 
WAVTAB 
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LIV. PROCEDURES 
we cael 


fee OBTAINING STEADY STATE 

Following program modification to reflect the geometry 
of the XR-3, program runs were made to obtain steady state 
conditions. This was accomplished by comparing steady state 
program outputs of thrust and pitch angle with measured 
@eweowootained from the XR-3 testcraft. The testcraft mea- 
Sured data consisted of steady state thrust and pitch angle 
measurements made at several velocities between 10 and 25 
knots in calm water. These measured data formed the basis 
upon which program modifications were to be made. 

It was first necessary to change the thrust versus 
velocity curve produced by the program so that it coincided 
with the measured data.- This was accomplished by adjusting 
the bubble frictional drag, ecoctf&leclent to bring the lower 
portion of the curves into agreement. A Similar agreement 
was obtained at the high end of the curve by adjusting the 
@emietrictional drag. This drag is a direct result of craft 
fedget, and 1s therefore not easily obtained. Adjusting draft 
requires that both plenum pressure and fan flow be maintained 
in such a manner that the appropriate lift is attained 
without exceeding the horse power limitation of the fan 
engines. Successive runs using the constant speed, trim 
option of the program were made at various drafts until the 


appropriate values of thrust for a velocity was obtained. 
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Pe iieeecuncmbniatmerarkt which produces the required 
eee, hence, the proper thrust, bubble pressure and stern 
Peale leakage were adjusted to bring the forces in the Z 
Seeeection into balance. To accomplish this the plenum 
pressure that produced the required lift was found and then 
the stern seal leakage was adjusted to provide the proper 
fan flow which maintained fan engine horsepower within its 
Mmits. A small portion of the total seal leakage was 
assumed to exist at the joints between the bow seal and 
Sidewalls, the remainder being under the stern seal. 

No attempt has been made to provide agreement between 
the craft transient behavior and that of the program. How- 
ever, the dataherein presented appears to show that there 
is a fair amount of agreement in this area. The transient 
behavior of the craft and the simulation program provide 


a very fertile area for further investigation. 


B. DIFFERENCES ENCOUNTERED IN THE MODEL 
meeceal Centiguration 
For reasons previously stated, the mathematical 

model of the 100B stern seal was selected to represent both 
seals on the XR-3. The forces ~ealculated by the simulation 
program employing this model appear to be reasonable. That 
is, they are of the same order of magnitude as those scaled 
from the 100B. Leakage coefficients were not changed since 
the type of orifices employed on the XR-3 were considered 


membe Similar to the LOOB. 
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2. Plenum Dimensions 

LO@pe mas wei at lNoLementing craft plenum geo- 
metry actual craft measurements were taken and used as input 
data to the simulation program. However, it was found that 
these dimensions did not produce appropriate XR-3 force 
and moment relationships. Through investigation of XR-3 
model data and subsequent program utilization the arbitrarily 
chosen nominal values of 10 and 20 feet, plenum width and 
length respectively, provided the desired results. Lift due 
to plenum pressure was highly sensitive to these plenum di- 
mensions. Therefore, when more accurate plenum pressure 


data becomes available these dimensions may require some 


ferinement. 


See eNSOITIVITY TO CHANGES TO INPUT DATA 

It was found that the XR-3 is extremely sensitive to 
the initial Be niteicns Pict Cutie eIipuL Gata Set.  omalt 
disagreement between the required initial conditions and 
ese Input, produce violent transients that can cause the 
program to stop execution or to give erroneous results. 

imerOousticmmost Critical input parameters are craft 
Graft and plenum pressure. This is true because the plenum 
pressure and draft are closely interrelated and cannot be 
separated for analysis. The plenum, having a large surface 
area, exerts the greatest single force on craft dynamics. 
Hence, a very small change in plenum pressure operating on 


this large surface area, 200 square feet, produces a significant 
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Saemge in the Z direction force. Thus, imbalance of either 
Graft or pressure, or draft and pressure, cause transients 
LOO great for the simulation to arcept. These results also 
Suggest that rapid maneuvers or inputs approaching step- 
functions will also produce transients too great to be ade- 
miacely computed. 

It is recommended that a investigation be conducted 
to determine methods of improving the integration procedure 


to allow very rapid maneuvers. 


De POSSIBLE ERRORS IN THE PROGRAM 

Many coefficients herein described were assumed to be 
correct. However, this may not be the case. Those coefr 
ficients that are calculated by their definitions such as 
Reynolds number and Froude number are most probable correct. 
Meee Caxage Orifice cocfficients, skin friction coef- 
ficients and the bubble drag coefficients have not been 
verified, though they seem to give adequate results. 

Other coefficients unknown to be exact are the moments 
of inertia. The simulation program provides a method for 
Seewrately calculating the moments of inertia, I T if 


and I This option was not employed because accurate 


a7 


weight distribution data was not available. The moments of 


Mmertia I WY felix I, Wenge GQoeeailnecerrom ceterences aes- 


Ko 
Beoing the craft prior to many major modifications. Not 
Pmiem react Of which was the installation of the auxiliary 


power supply and the replacement of bow and stern seals. 
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These modifications added a considerable percentage to the 
total craft weight. These weight additions were concentrated 
in several locations. Hence, the scaled values of Ty. IY and 
I, used in the simulation cannot be relied upon as being 
mecurate. 

No available data could be found concerning the product 
moment of inertia Tye This moment of inertia 1s required 
when attempting to execute maneuvers involving turns or off 
centerline craft motions. Therefore, an estimated value 
using Size scaling was employed and rendered reasonable craft 
motions under the above conditions. 

To determine the correct moments of inertia, it is re- 
commended that a thorougn craft mass distribution analysis 
be conducted and the resulting parameters be input to the 
program omens the rer mentioned above. 

A third area open to further improvement is that of 
refining the dimensions input to the sidewall map of the 
side geometry subroutine. The present representation for 
the XR-3 sidewalls loses much computational accuracy Gaue Lo 
the Laneee vertical stepsize presently used. 

A remapping of the sidewall geometry is required to 
more accurately describe the curvature in the lower regions 
of the sidewalls. This would result in the subroutine pro- 
viding more refined values of sidewall volume and surface 


Peeaetor the Various XR-3 drafts and attitudes which are 


Miportant for the calculation of sidewall.forces and moments. 
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Notwithstanding the above areas in which error exists 
teers felt that the craft data and figures tend to indicate 
that the simulation program is a valid representation of 


XR-3 motion. 
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Vio nkocmanM Viki Fr TCATION 


foe LNSTRUMENTATION AND RECORDED DATA 
Instrumentation as presently installed in the XR-3 
eludes the following: 
1. Two strain gauges mounted directly on the engine 
mounts measure the thrust of the propulsion 
SS cen. 
2. A rate Gyro mounted near C.G. is used to measure 
pitch and roll angles. 
3. A turbine type flowmeter which extends forward 
in undisturbed water measures velocity. 
Additional instrumentation presently operational on 
the XR-3 meaSures yaw, pressure in the plenum and seals, 
and rudder angles. The. output from these sensing devices 
are recorded on a 14 track Pemco Model 120B Magnetic tape 
recorder. Hard copies of the measurement data are obtained 
uSing a Hewlitt Packard two channel chart recorder. No 
filtering other than that natural filtering caused by the 
frequency response of the recorders was used with the re- 
sult that the data was quite noisy. The noise was consi- 
dered to be white and gaussian, therefore the mean value of 
the resultant curve was taken as the best estimate of the 
Variable under cosideration. 
Shown in Table I are the data collected from test runs 
Mage by the XR-3 at Lake San Antonio, California, over a 


Sne-day period in calm water and no wind conditions. Each 
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mepresenes a Steady state condition. 


uy 











Patoile £. 


Velocity in kts. 


24.5 


Ee. 0 


Jiabao bichon ot elerbh eters 


405 


A25 


S)40) 


475 


453 
490 
400 
465 
407 
Bo 
415 
425 
440 
455 
460 
455 
480 
496 
SLC 


SOU 
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Test Craft Measured Data 


Pitch angle in degrees 


NOt available 


N@twuavallLeable 


cps rs, 


Not available 


Not available 


Not available 
Not available 
Not available 


N@travaalaple 


Not available 
Not available 
Not available 
Not available 


Not available 


Neieeavadl lable 








B. L&M PROGRAM DEVELOPED DATA 
aly Seach Sem usein Conditions 
Steady state trim conditions for the XR-3 using 
the L&M program were difficult to obtain. It was observed 
that the program initial conditions must be very close to 
the steady state values or the imbalance in forces generated 
within the program iS so great that one of two cases occurs. 
Gace using Variable Time Step (R-K—M) 

The force imbalance causes such a radical transient 
mimic, In order not to lose curve information, the inte- 
gration step size is reduced until the minimum allowed step 
Size is reached and the program execution stops. 

G@oewe using Fixed TL) eneises (ON) 

The danger in this method is that the step size may be 
too great (note: in Sone cases, .01 seconds is too large) 
and curve Miteunacron 1s lost. When information is lost 
the program calculations begin generating erratic number 
values thuS causing an underflow or overflow condition to 
exist within the program which stops program execution. Thus, 
care must be exercised in selecting a stepsize to be used. 
The results shown in Table II were obtained using the fixed 
step size option with a step size of .001. The constant speed 
S@emevatiable thrust options of Block 1 of the INCON Subroutine 
were selected. Since the force equations of the program are 
dependent upon speed, using these options removes the surge 
Equation from the computations and adjusts thrust to balance 


mie L£orces in the X direction. 
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Table II. L&M Program Initia] Conditions 


Speed Thrust Pitch Angle Draft 
(in knots) (in pounds) (in degrees) (in inches) 
LLC) at SOS eee 6E16 
2.5 4 AS: sal 6.45 
ns, 0 BS) oe 94 62.55 
i 5 Bey © 6 655 
Zor 0 424.7 29 5208 
Bee} 464.7 lies: 52.18 
Zor 0 bez 28 09 Bio 
a SOG 8 e . 06 5a 3 
30.0 625.5 Bons: Bao 


Plenum pressure 24.86 psf 
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To balance the forces in the Z direction short time 
runs Of .01 seconds or less were used. For each run ad- 
justments were made to the initial draft and/or plenum 
fMeessure tO bring these forces into balance. Care was 
taken to avoid fan flows and pressures above that level 
requiring greater than 15 horsepower, which is the maximum 
horsepower available from the fan engines. Glo) en ealielslabeve fea gl 
this manner the initial trim conditions of draft, Dateeh 
angle, and thrust required were obtained for each of the 


speeds shown in Table II. Once these initial conditions 


were known, the step size was increased to .01 and it was 
momma that no curve information was being lost indicating 
that the transient was of a less violent nature. Long time 
runs uSing the variable time step size verified that these 
values eee the COmEeecte Mntiidie titii@scOnarttions. A Straight 
line segment was used to join the points so that linear in- 
terpolations may be used between the tabular values. 

To obtain good agreement between measured and simulated 
data required the manipulation of two quantities that pro- 
emce total ship drag; the bubble drag, FXPWAV, and the skin 
fetetion drag FX. Adjustment of the forces in the X direction 
was accomplished by adjusting draft to provide the proper 
Meag force at the upper end of the velocity curve where the 
frictional drag component exerts the ase Biiuience Overs 
Bera!’ Ship drag. The next step was to adjust the coefficient 
Seetnie bubble drag to provide the minimum point in the total 


drag curve at the velocity indicated by Table II, about 17.0 
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maees, and to also provide the required drag at the low end 
Memene Grag curve where the bubble drag exerts the greatest 
influence. 

To obtain agreement between measured and simulated 
pitch angle required some programming changes. Reference 11 
indicates that a pressure gradient exists in the plenum 
chamber. This gradient is such that as the craft increases 
in speed the center of pressure moves aft. This causes the 
craft to have a pitch up attitude at slow speeds and a gra- 
dually smaller positive pitch as speed increases. Using the 
experimentally measured data shown in Table I, a set of points 
representing the craft center of pressure position as a 
function of speed was calculated which would produce the 
appropriate craft pitch angle for a given speed. A linear 
function used to approximate these points was then included 
in the simulation program. 

Figures 14, 15 and 16 show plotted points of the two 
tables of:data. Triangles represent experimentally deter- 
mined data and squares Pew pregraMm Output. Clearly the pro- 
gram Simulates quite well the thrust and pitch angle curves. 
mecurate data on draft was not available from XR-3 tests 
due to a defective height sensing device. However, it 
memecelt that the draft data when available will closely 
fit the curve shown. 

ue eee aneeton OL Teed Program Data 

Figures 17, 18 and 19 demostrates the use of the 


feerapble thrust Option. Thrust is input to Subroutine PROP 
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Figure 14 
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Figure 16 
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meebizing the Block 16 variable thrust option of Subroutine 
miecON. This thrust is applied equally to both screws and 
is mapped as follows. Thrust for 20 knots is initially 
applied and maintained for 10 seconds. Between time 10 
seconds and time 15 seconds the craft experiences a ramp 
meme tion Of thrust building up to the tabular value for 

25 Knots. This thrust is maintained for the remainder of 
the run. 

Figure 18 shows that the pitch angle follows the ex- 
pected curve, in that initially pitch up is experienced 
as additional thrust is applied. This is followed by gra- 
Gually dampening oscillations toward the value of pitch 
meme expected at 25 knots. 

Figures 17 and 19 show that speed through the water 
increases and that draft decreases, aS expected. Addi- 
tionally, examination of Block 16 of Subroutine INCON makes 
it clear-that the engines may be unevenly thrusted rather 
than equally thrusted as shown. 

Figures 20, 21, 22 and 23 show the behavior of the 
eran for a turning maneuver. The rudder input vs. time 
is shown in Figure 20. Corresponding roll angle vs. time, 
Pitch angle vs. time, and the ship horizontal plane motion 
are in Figures 21, 22 and 23 respectively. 

Pemuomincehesting to note that the roll and pitch 
aeamping coefficients appear to be different during the time 
the rudder is off centerline than when the rudder is re- 


turned to centerline. However, these differences are to 
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berexpected due to the nonlinearity of the system. The 


Soccillation appears to be about 0.5 H This oscillation 


Z° 
was not felt by the authors during craft test runs. However, 
it may not be sensed due to its high frequency and small 


fmemetude. Unfortunately at this writing no accurate turn 


data iS available. 
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fie DLoCUSo LON AND EVALUATION OF SIMULATION RESULTS 


A. PioCUSSION OF PRODUCT MOMENT OF INERTIA; Tyg 
Peminitormatiton Mas beensavailable on the value of the 


product moment of inertia I fOr the %R=-3. Therefore, 


XL 
throughout the XR-3 program development the moment of inertia 
lio was set equal to zero. This decision had no effect on 
meate MOCIOnS for calm water Straight runs Since craft 

motion was restricted to translational movement in the Z 
mae. direction and rotation about the Y axis. However, 

when turns were attempted and waves off the bow were intro- 
duced, program instabilities were observed. These may be 
seen by referring to Figures 24, 25 and 26 which show that 
for small perterbations the craft reacts in a stable manner, 
but as the ruddér angle increases toward 20 degrees the craft 
becomes unstable even though waves have not yet been intro- 
duced. A weight scaled value of Iya was obtained from the 


100B input data set and was later used as a representative 


merue for the XR-3. 


B. 5 DEGREE TURNS WITH AND WITHOUT WAVES 

nan Seeenoe Gemverlty —COMNnGlel1Ons reported in Reference 
Oetor the 100B model, a 20 knot run with a 5 degree rudder 
turn with and without waves was introduced. The wavelength 
and wave height were scaled from Reference 8 according to 
relative craft size. Hence, a wavelength of 72 feet and a 
wave height of 0.16 feet were uSed with the XR-3 model. Wave 


amplitude and frequency of encounter are shown in Figure 29. 
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It is interesting to note that this same frequency can be 
seen in roll angle and pitch angle as shown in Figures 27 
and 28. it is observed in Figure 27 that the amplitude of 
roll increases almost linearly as the wave front moves away 


from the ship's head, as would be expected for beam seas. 


meer uURNS WITH 5, 15 AND 20 DEGREE RUDDER 

Mmieites 3), 32, and 33 show the roll, pitch, and yaw 
motions produced by 5, 15, and 20 degree rudder turns at 
20 knots in calm water with the product moment of inertia 
Ivy inserted. In all cases, damping and system Stability 
ere observed. Similar turns made with Ivy equal to zero 
revealed the presence of roll and pitch instability for 
large rudder angles. In both cases where 5 degree rudder 
turns were made the variations in the roll and pitch angles 


Smeaned small and decayed to zero within 15 seconds or 


mieout 5 cycles. 


D. COMPARISON OF SIMULATION AND MEASURED DATA 

At the time of this writing meaSured pitch and roll 
data was not in a form suitable for definitive comparisons. 
However, recorded data of testcraft motion during turn ma- 
neuvers as shown in Figure 34 tend to validate the presence 
Smecamped roll and pitch oscillations which are also pre- 
sent in the simulated turns as ser iMmiacures si, 32 and 33. 

Measured and recorded plenum pressures tend to cor- 
respond with those calculated by the simulation program. 


This pressure appears to vary between 24-25 pounds per square 
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Figure 30 
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Pigute 32 
FIVE, FIFTEEN AND TWENTY DEGREE TURNS: 
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Foot throughout all the various monitored craft motions, 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


The conversion of the L&M program from the 100B model 
to the smaller XR-3 model involved considerably more time 
and effort than was anticipated at the onset. A great deal 
of time waS required for familiarization with the program, 
Por PLrOViding a more useful output format and for adding 
plotting routines. Considerable effort was required to 
Mieeain the straight run, calm water trim conditions for 
various speeds. Plenum center of pressure movement with 
velocity was one of the more difficult problems encountered 
mmootaining the proper initial trim eocee con oan ileal} 
recommended that future investigators write an initial 
condition program to assist in reducing the time and effort 
reguired to arrive at the proper trim conditions. 

Seousagreemena between measured and computed data for 
calm water straight runs was possible by the adjustment of 
the hustle Stag seOctilcienes, Gdratt, bubble pressure and 
stern seal leakage. The meaSurement data for turning ma- 
neuvers could not be used for a check on the simulation in 
meee tOrm In which i€ iS presently recorded. It is recom- 
mended that the recorded data be processed digitally and 
Suitably reduced to a form useful for comparison with the 
L&M program output. Qualitative comparison of the XR-3 turn 
maneuver data however, haS indicated that pitch and roll 
motion have a lightly damped response quite similar to the 


L&M program output. Much work remains to be done in checking 


EOS 





the accuracy of the computed transient response. 

Additional qualitative checks were made on the con- 
puted XR-3 turn maneuvers by comparison with the 100B 
model response. The results indicated similar transient 
Saeracteristics in which damping decreases as pitch or 
meml angle increase. 

The simulation study of the XR-3 has demostrated the 
importance of the product moment of inertia Log in main- 
taining stability in roll and Puech fOr ridder angles 
meeacer than 5 degrees. Future investigators should attempt 
to obtain an accurate value for Togs This can be done by 


utilizing the discrete mass distribution SOGrome Liltiic 


INCON Subroutine. 
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FE RODUCTION 


The SES Motions and Loads Program is a Fortran IV pro- 
gram capable of calculating the time histories of the motions 
and loads of a surface effect ship operating in the on bubble 
mode in a seaway in six degrees of freedom. It is programmed 
ina "modular" fashion,i.e. each major element of the craft 
has its own subprogram. The input is in the form of basic 
data, i.e. sidewall offsets, seal PVeONerry ean map. data 
points, etc. and enables the user to make a wide range of 
parametric changes to the particular craft being modeled 
by changing only the input data. The program has the ca- 
Pability of printing and/or plotting a wide Bane of var- 
jables as a function of time; such as forces and moments 
from each element, total craft forces, moments and ac- 
celerations, leakage gaps and flow Bacco, tan tlow rates 
and power aaGieee eee Uaewecurrecnt program has been 
checked out and run on an IBM 360 CoMpunzer (using version 
Release 20.6.) at the W. R. Church Computer Facility, NPS, 
Monterey, Ca. 

This report is intended to be primarily a users manual 
BOr the program as it exists. The details of how one would 
rewrite various modules to reflect changes in the craft de- 


Sign are discussed extensively in Appendix A-~A, 
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PROGRAM ORGANIZATION 


Within a given run the program proceeds as follows. 
First the main program calls the input and initialization 
Subroutine INCON. The program then proceeds to calculate 
the motion time histories by calling the various modular 
Subroutines and the integration routine for aS many time 
steps as is necessary to complete the run. During the 
course of the calculation the various subroutines will 
output those variables selected by the print option switches 
Memeene appropriate input print interval. In addition, 
data may be written on Scratch files to be used for sum- 
Marizing the output and doing bending moments and shear 
Calculations. At the present time the Shears and Moments 
Pw@eroucine has been removed, however, it may be reintroduced 
at any time with a minimum of changes. 

When the run is completed, the program, depending on 
the options chosen, may then print and/or plot the output 
summary and proceed to do the bending moment and shear cal- 
Sulations. After ee steps are completed the main pro- 
gram returns to INCON to read data for the next case. 

The above procedure 1S repeated until all cases are 


finished. , 


ie 





PROGRAM REQUIREMENTS 


Asecun on an LBM S60 computer under OS/360MVT version 


(Release 20.6), the program requirements are as follows; 





LOA ree ORE TIME 
Program length 201,456, NA 
Core Size required for execution 250,000, 4 SEE 
BELOW 


The total core required for a given step includes the 
program, system routines, data storage, etc. The time re- 
quired to run a single case depends on the nature of the 
Case and the length of time to be simulated. The ratio of 
required execution time to real time has been found to be 
as low as 1/3 for some calm water cases. Higher ratios 
meen where the transient curves are rather sharp. 

The program currently uses the following FORTRAN Data 


Set Reference Numbers (Unit Numbers): 


Unit Number Usage 


i Pematenerile For COLFIL plotting package. 
2 . Seroavecneeriic efor Lateral plane output 
summary 
5 Input 
6 Oulpmt 
GAS Sidewall Integral Table data (see Appendix A-C) 


ae 





IBMF Intermediate file containing time histories 
required by bending moment calculations. The 
unit number IBMF is read in using the 01501 


input card (see discussion of Input) 
Bieo AN OULPUT 


The input for the program is handled by subroutine INCON 
and is discussed there. The output is controlled by the 


individual subprograms and is discussed separately there. 
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DESCRIPTION OF SUBPROGRAMS 


eS) 








Main Program 


DiscusSion - The main program contains the logic for inter- 
mennecting the various subroutines. In addition, it com- 
pares the running value of time with the finish time, 
Calculates the next value of time for printing after each 
print time, does some output, etc. The main program also 
Peneains the logic for the integration of the trajectory 


(x, y,WV) by the trapazoidal rule. 
Input - None 


Output - The output from the main program is controlled by 
meee print Option switch ITRAJ. If the value of this switch 
is 1, the program will print the current values of the time 
meiamethe Craft velocities and displacements (angular and 

translational) as shown in the sample output at the end of 


the report. 
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Subroutine AEROD 


Discussion —- Subroutine AEROD calculates the aerodynamic 
forces and moments on the craft. The subroutine currently 
uses curve fits of the aerodynamic coefficients vs sideslip 
angle. Appendix AA discusses how the program may be changed 


to input tabular data for the aerodynamic coefficients. 
Usage - CALL AEROD 


Mmeput — The input data for AEROD is controlled by block 10 


miesubroutine INCON and is discussed there. 


Sutput —- The output from AEROD is controlled by print option 
@e'echn LABROD. If this switch is 1, the program will print 
the total forces and moments calculated as shown in the 


Semple Output at the end of the report. 


eel 





Subroutine BOWSL 


Hascussion —- Subroutine BOWSL calcualtes the forces and 
moments due to the bow seal and the leakage flow rate 
associated with any gaps which open up under the various 


Seal stations. 
Usage ~- CALL BOWSL 


Input - The input data for BOWSL is controlled by blocks 


feo, and 1/7 of subroutine INCON and is discussed there. 


Seeput - The output from BOWSL is controlled by print option 


Switch IBOWSL. If this switch is 1, the program will print 


the gaps and wetted lengths at the various stations on the 


bowseal as well as the total forces and moments due to the 


bowseal. A Sample print point iS given in the sample out- 


put at the end of the Ge port. 
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embroutine COLFIL 


Discussion - Subroutine COLFIL is used to format the pro- 
Gram output. It can provide up to 16 variables listed in 
two tabular summaries of 8 variables each. The subroutine 


has a plotting package that permits use of either the print 
Brot Or Ccalcomp (continuous line) option. Multiple curves 
may be plotted and up to 10 graphs may be drawn on each run. 
Print switches control the printing by each subroutine of 
its forces and moments. This option is very useful in de- 
bugging the program. Normally it provides too much output 
and is too cumbersome. 


Usage - CALL COLFIL 


Mmeue — The input is controlled by the control card with 
ime control tae 00104 and is discussed in the description 
of Subroutine INCON in the block 1 control card section. 
ParcLonal input required for the plotting package is con- 


Meemeodein the Users Manual for the COLFIL routine. 


Smeput — For output format consult the COLFIL plotting 


package users manual. 
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Subroutine DMINV 


Discussion -— Subroutine DMINV is a matrix inversion package 
used for inverting the mass matrix. The method used is the 


standard Gauss-Jordan technique. The determinant of the 
[errs 1S also calculated. If the determinant is zero the 
matrix 1S singular and the program will stop (see Appendix 


A-B Abnormal Termination Messages). 
Usage - CALL DMINV (A, N, D) 


Description of Parameters 
A = The input matrix, destroyed in computation and 
replaced by resultant inverse 
N = The order of matrix (maximum = 6) 


mee rie resultant aeterminant of the matrix 
Input - None 


Sireout — None 
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Subroutine FAN 


Miecussion - Subroutine FAN is used to calculate the inflow 
to the plenum, stern seal, and bowseal from the fans for 
a given pressure differential across the fans. The sub- 
routine uses tabular data and allows for changes in fan 
Speed as well as the number of fans. Fan speed changes are 
currently handled by using the fan scaling laws rather than 


reading data as a function of fan speed. 
weage -— CALL FAN 


Sent -— The input for subroutine FAN is controlled by block 


moor SUbroutine INCON and is discussed there. 


Satput — The output from FAN is controlled by print option 
Switch IVEL. Ihe Wy elolaic Fea ch is 1, the program will print 
the flow rates and pressure differentials across the bow 
seal, cushion and stern seal fans, as well as the total fan 
horsepower - required as shown in the sample output at the 


emer ot the report. 
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Function FGL 


Session - The purpose of FGl is to perform the task of 
evaluating a tabular function y = f(x). The procedure used 
is linear interpolation between data points. If the input 
value of the independent variable does not fall in the range 
of its input tabular values the appropriate table extreme 
is held. The function y must be single valued and the ta- 
bular values for the independent variable must be stored in 
ascending order i.e. XTAB(J). XTAB(J+1). The increment be- 
tween tabular values need not be constant. 
Meage —- Y = FGl (X, N, XTAB, YTAB, IMEM) 
Description of Parameters 

Y = The output value of the tabular function 

X = The value of the independent variable to be used 

teomevyalvate the function y = £(x) 


N = The number of elements in the XTAB and YTAB arrays 


~ 
i 
a 
w 
II 


The name of the array that contains the values of 

the independent variable in ascending order 

YTAB = The name of the array that contains the values of 
the dependent variable corresponding to the values 
of the independent variable stored in XTAB 

IMEM = Dummy variable used to store location of last entry 

imeem taple, Btoube Used as initial point of search at 


erce went a VY 


Input - None 
/ 


Sutput —- None 
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Subroutine FORIT 


Discussion - Subroutine FORIT is used to obtain a Fourier 
analysis of a periodically tabulated function. It computes 
the coefficients of the desired number of terms of the Fourier 
Series £(x) = a(0) +B (alis)cos Poy oii KX), LO approx= 
imate the given set of periodically tabulated values of a 


mmctilon. 
Seraege) - CALL FORIT (FNT, N, M, A, B, IER) 


Description of Parameters 
Pmitis== Vector of tabulated function values of length 2N+l 
N = The interval such that 2N+l points are taken over 
PiemenecervalmO, 2). The spacing is thus 2Tl/(2N+1) . 
N must be greater than or equal to M 
M = Maximum order of harmonics to be fitted. M must 
be greater than or equal to zero 
A = Resultant vector of Fourier cosine coefficients 
B = Resultant vector of Fourier Sine coefficients 
IER = Resultant error code where: 
ITER=0 No error 
IER=1 N less than M 


IER=2 M less than 0 
iput —- None 


Output - None 
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Subroutine INCON 


iscussion oo soummmestNCON Contains the logic for the 
meaoing Of all input data, the initialization of variables 
and for the initiation of new cases. The subroutine is pro- 
grammed so that the input data is non-structured i.e. there 
is no set order of input cards from one read statement to 
another. However, a given read statement may have more than 
one input card. In runs with multiple cases, only those 
input parameters which are to be changed for the new case 


meea be read in. 
Beage — CALL INCON 


aiewe =~ Lhe input logic is subdivided into twenty-two “blocks. ' 
Meena olock being for a particular subroutine or function e.g. 
stern seal, sidewall, etc. Switching to the proper input 
i9lock is done by means Of a control tag on the input cards. 
Smee Control has been transferred to a particular block ad- 
@eeLonal input may be read by specific format statements if 
necessary, however, = CaeGmviteneecontains the control tag 
may also be used to input up to seven variables for the 
selected block. The format of the input card used for swi- 


mening is as follows: 


Column Format eee 
Ss a5 SOnerel tag 
4-5 T2 option tag 
OS Ies: 7F10 mmememvactables (if any) 
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iiewecemeroletagiiscea FLO transfer control to the appro- 
priate block is of the form OMN where M and N are integer 
PemstantsS, 1.e., if the control tag is 015 the program will 
MO tO Statement 1500. 

The option tag is loaded into a variable called IOPT 
and allows different branches of logic within a given input 
mock. Lf no branches are used within the input block the 
option tag may be zero. 

The input variables are stored in an array called TEMP 
and may be set to a given parameter by uSing an arithmetic 
replacement statement in the appropriate block (see various 
blocks in attached listing for examples). 


Example: When read, a card of the form 


Column So a ass ie Sie) 80 
¢ b dy & 4 
O17 02 B25 blank iene) blank 


will transfer control to statement 1700. Within statement 
1700 the logic might look as follows 
700 CONTINUE 
Comoe MOrm 720), OPT 
iv lt eA = TEMP (1) 
B = TEMP (2) 
GO TO 10 


720 Ce EMP (1) 


D TEMP (2) 


is TEMP (3) 


COO LO 


aS) 


In the above example the option switch IOPT is 02 and 
will cause the program to go to statement 1720 which will 
mem set C=3.5, D=0.0, E=1.0 and then will return to statement 
10 which is the read statement for the input cards. 

The various input blocks as currently programmed are as 
follows: 
Block Number Function 

i. Executive input (time steps, print switches, 


option switches, etc.) 


2 Mass and inertial input 
3 Craft geometric input 
4 Sidewall input 
5 Stern seal input 
6 Bow seal input 
5 Plenum input 
8 | Propulsion geometry input 
SS Rudder input 
Ie, | Aerodynamic sl ayenee 
lb Wave input 
12 Peer imcond li lons 
ys biecesOn case, Indicator 
14 Ener lydreator 
1S Not used 
16 Pieist  wuddem mOE1On Lnput 
lg Bow and stern seal pressure input 
18 ae Cease 
ILS Fan map input 
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Card Columns Bormart acy, 


u ka T DOdOiS= control tag and option 
Ea ILS F Start time (generally zero), sec. 
16-25 F Hats tame. Sec. 
NDS S. F Initial integration time step, sec. 
36-45 F Print interval, sec. 
46-55 F Print start time, sec. 


Option tag = 2 Palieswitones  (1f=l print if=0 no print) 


Card Columns Mousse BME ry 
HL ah 5 af HOm02. — control tag and option 
2 5 a HWACeri, souint switch’ tor ,Lateral 


acceleration and integration time 


SHEEN ®. 

10: ioe. IVEL, print switch for subroutine FAN 

15 1) Deed eorink Switch for main program 

20 I Lover switch for subroutine 
SIDEWL 

2m i IBOWSL, print switch for subroutine 
BOWSL 

30 a iow ist esogiinien switch tor subroutine 
STNSL 

35 1p IWAVES, print switch for subroutine 
WAVES 

40 I ike print switch for subroutine 
RUDDER 
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Motion tag = 2 (Continued) 


card Columns Format JB) geen ay's 
2 45 I PeROEy Dicintk switch for sub- 
Meont. ) 


routine PROP 
50 al PAEBROD print Switch for ‘sub- 


routine AEROD 


55 it iow oaiMt switch tor Subroutine 
RHS 
Option tag = 3 Integrator Tolerances 
Card Columns Format Entry 
it Les ip OOMOe  recOnEerol “cag and Option 
2. 1-2 1, Number of integrators to be used 


(maximum of 15) 


3-4 LZ Somuroletag fOr Jitegqration pro= 
cedure 
ed. . « oO 8F1O0 Pieeaq@aton Crror tolerance, one 


per integrator 
Note: The usual procedure for selecting integrator tolerances 
is to choose values a the order of 107° to Looe However, 
m@eieweach new craft the following procedure is used; make a 
run with a given set of tolerances; halve the tolerances and 
Zone Lt the results do not change drastically, try doubling 
the original tolerances and running. As changes are noted in 
Peewee eGiatOor Outputs (1.e. craft velocities, etc.), hold the 


tolerances for those integrators fixed and vary the others 


until the largest tolerances possible are obtained. Once 
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chosen for a given craft the tolerances Mawely Nave to. be 
changed. A control tag (JQQ in the program) is used to se- 
Meet integration procedure. A blank provides the original 
R-K., procedure. A"one" provides a fixed step size runga- 
Kutta procedure that could be useful in selecting integrator 
tolerances. A "two" provides runga~kutta variable step size, 
but does not print out the step varying procedure. Org choles 


minimum stepsize is Oe out - 


Beevon tag = 4 Tabular Summary of Output 


ard Columns POrivat Bite 
ie 1-5 il 00104, control tag and option 
2 5 if Print switch for summary one 


aoe ey Ol Oe int 
10 io Print switch for summary two 


iL = jeheatone Atel eaves opens’ ol 


Option tag = 5 PEOGEamM Option Switches 
eard Columns Format Entry 
1 1-5 | I 00105, control tag and option 
6-15 F Switch for lateral plane motions 
Poul lateral plane motions 
O= six dégrees of freedom 
16-25 F wiecitOr Constant speed 

1 = surge equation not used 
0 = surge velocity allowed to vary 
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eard Columns Format Entry 


26-35 F Svairch for tram 
1 = thrust varied to give proper 
Se-Ccdee In Enis Opticon is se-— 


PeCKCGysseNeh. there MUSE. De anpute 
in Block 16 an estimate of thrust 
and a 0.0 for the number of data 
ole mbie hae 

J°—-SEhrust as input 


Block 2 - Mass and Inertial Input 


Sption tag = 1 Summary Mass Properties Card 
eard Columns Format Hey, 
L Los dl O0Z Ota seOntre ly tag 
OS Irs F- Total craft weight, pounds 
16-25 F EOnGgmeudinal Center Ob Gravity, 


feet forward of transom 
26-35 F Vertical center of gravity 
feet above baseline (keel) 


Bio-45 2 I 


xx Mass moment of inertia 
pout cds = Slug — ft. 
46-55 F Toy! mass moment of inertia 
Poeottnyeds1S, Slug — Bee 
56-65 F Loot mass moment of inertia 
ap9oue Z axis, Slug - ft.? 
66-75 FE. Le; MaAssmmement Of inertia 
deome xz ax1s, Slug - ft. 
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Note: If this option is used, shears and moments will not 


Mewcalculated correctly. 


Metion tag = 2 Bisciere MaSs Distribution Cards 
Card Columns Format Entry 
1 LS Ai W202 = control tag 
feo, 4... 1 -10 F Discrete weight, pounds 
20 F Longitudinal c.g., feet forward 


Or Eransem 


bo 0 F Transverse c.g., feet to starboard 
31-40 F Vertical c.g., feet above baseline 
(keel) 
Last ie) F iO CON tO. lag 


eanrd 

The mass distribution is specified by discrete weights at 
various iSeartseone with ae card used for each weight. Sym- 
metry about the longitudinal center-plane is assumed, so that 
only weights on the eeeeeeaed Side, with atransverse location 
Meedeer than zero should be specified. Thus, the total weight 
Sieene input cards is one half the total craft weight. The 
last Ed is used to signal the end of the discrete weights. 


The maximum number of discrete weights is 201. 


Block 3 Crake Geomerric Tnput 


No option tags 


Card Columns Format Enea 
ii L=5) i Op 00m econerol tag 
6=15 EF NUM@pemnOme Stations on port Sidewall 
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Card Columns Senaiialias inet 


16-25 F Number of stations on starboard 
Sidewall 
26-35 F Number of stations on bow seal 
36-45 F Number of stations on stern seal 
46-55 F hewwlmciakt ength, Lt. 
Note: The maximum number of stations is 11 for each element. 


The program will automatically subdivide the element into 


Semesections (N is the number of stations) and calculate the 


xX and y coordinates of each station. 


Block 4 Sidewall Input 


Option tag = 1 Sidewall Input 


Card Columns Format Pittry 
1 lan, T 00401 - control tag and option 
eas F y distance from centerline to 


Sioewal ls. ££. 


i G20 | F Average wetted length of side- 
| Welded tee 
26-35 F neaivagie otlrtce COoecLLicient of 
Sidewall 
BG-d5 18h Gtece-itlow drag coerficirent of 
Sidewall 
46-55 F Average beam of sidewall, ft. 
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Meee: Under this option tag the program will also read the 


Sidewall integral table file stored on Unit 10. 


Brock 5 Stern Seal Input 


No option tags 


fara Columns Format yatta ay 
i ees. a HOS00F = Control tag 
1 6-15 F MecOOradinace Of Seal hinge, £t. 


LOpwacd. OL cCransom 


16-25 F Zecoordinate Of Seal hinge, £t. 


above keel 


ID ca 3) F Base leakage area, sq. ft. 
Ihara ae Seal leakage orifice coefficient 
46-55 F Angle between leading edge of seal 


anewemarme Vert:cal, deg. 
56-65 F Pressure dirtterential between 


seal bag and bubble, psf 


66-15 lah Length of leading edge of seal, ft. 


HS? 





Bock 6 Bovweoeal Input 


No option tags 


Card Columns Format Entry 
i AS al COGO0F-Seontrol tag 
6-15 F xX coordinate of seal hinge, ft. 


forward of transom 

Ie S) F Seal leakage orifice coefficient 

EASA E IS, 1h Pressure differential between seal 
bad sand bubble | pst 

36-45 F Z coordinate of the seal hinge 
above keel, feet 

46-55 F Angle between leading edge of 
seal and craft vertical, deg. 

56-65 | F lemetioog: the leaging edge of the 
seal, feet 

66-75 EF Base leakage area, ft. 


Breck 7 Plenum Input 


Set10on tag = il Plenum Geometry 


eS 3 


te 





Card Columns ong ene Entry 





i. aS I iy wie cCOntrol tag and Option 

il 6-15 F Plenum length at water surface, ft. 
16-25 F Plenum width at water surface, ft. 
26-35 = Not used at present 
36-45 F Plenum width at deck, ft. 
46-55 ip Peewume lenctin at deck; ft. 
56-65 F MecOoordinate:. Of center of pressure, 

Lew tOorward Of “transom 

66-75 F Plenum average height 


weeLion tag = 2 Critical Froude Number 


Card Columns Format Rimes 
1 ssi) Ie (vee control tagrand Option 
ba 1s F Froude number corresponding to 
hump speed 


meen S © Propulsion Input 


OoELOn tag = 1 Basic Data 


Card Columns evemMa t ghee 
IL eA F Hic ocontrol Eag and option 
CIS F Seeceoeorcimate Of propeller center, 


tEomTara Of LranSon 

16-25 le wecwctance from centerline to 
paome Willer Genter, It. 

Zoo F Zecoocemmace Of propeller, ft. 


above keel 
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Eeeock 9 Rudder Input 


Mpoe1on tag = 1 Rudder Geometry 


Card Columns Format Ha wis 
1 as) AE VUE = COnerol tag and Option 
6=15 F x coordinate of centroid of rudder, 


tes 2Orware OF transom 

16-25 F Vecistance from centerline to rudder 
Generold, ft: 

26-35 F ZeCe@Gnolnare GEyCentcrola Of rudder, 


ft. above keel 


36-45 F Rudder span, ft. 

46-55 F Rudder aspect ratio 

BOSS. F Ru@@er area, Sq. L£tx 

SG— 75 _ Average Thickness ratio of rudder 
| section 


mitock 10 Aerodynamic Input 


No option tags 


Card Columns Format meta, 
1 - Ihe 5 1 COO — Concrol tag 
6-15 F Reference length, ft. 
6-25 F Reference width, ft. 
Note; The reference length and width are those used to 


nondemensionalize the wind tunnel data used in the program. 


Block ll Wee Sa lmpue 
Option tag = 1 Amplitude vs. frequency 
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eard Columns 


u 5 


Seles) 


Om 5 


meee) 4... 1-10 


T1-20 


Option tag = 2 


FOrmac 


II 


Le 


E 


i 


hic ry 


Tee ecOntcrol tag and Option 


Number of wave components, max- 


imum Ot 10 


Initial wave heading relative to 


craft headings (180° = head seas), 


deg. 


Wave frequency, rad./sec. 


Wave amplitude, £ft. 


Amplitude vs. wave length 


Card Columns Format Entry 
1 i=-5 I Hie es —-cCONcrOl tag and option 
6-15 F Number of wave components, max- 
imum of 10 
16-25 F Initial wave heading relative to 
craft heading (180° =: head seas), 
Geq- 
meee 4... 1-10 FP Wave length, ft. 
ILS Z10 F Wave amplitude, ft. 
noreee The program constructs irregular seas by adding 2oO— 


gether a series of regular wave components with appropriate 


floret ouelon Gf amplituae and £recuency. To run a regular 


sea case, the number of wave components should be equal to 


Mmtimenounumber Gf waves 1S equal to zero (calm water) 


m@mwaadi tional data (cards 2, 


3, 


Block 12 initial Genditions 
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No option tags 


Card Columns FOrmat Entry 
i = 5 i OEZ0UF — control “tag 
6-15 F Initial speed, knots 
16-25 F Initial pitch angle (positive- 


nose up), deg. 
26-35 F Prieto orate Ene ‘Center .OL 
gravity, inches 
36-45 F Initial bubble gage presSure, psf 
mleck 13 End of Case 


Ne option tags 


Card Columns Format ie ite tay 

I Is i Os Ve Control tag 
Mete: This card is used to Signal the end of data fora 
Given case. The program Wihenow tnhitkalaze all required 


Meera bles and will start the calculation of the time history. 
Tf the program reads two consecutive 01300 cards with 
no data between, the job will be stopped rather than repeat 


the previous case. 


Block 14 Pncesor Run 


No option tags 


Card Columns Format fete Tey, 
i L=-5 iL Clavue— control tag 
Note: This card is used to signal the end of the run. When 


encountered the program will print out the message: 
COMPLETED ALL RUNS 


and will then stop. 
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mock 16 Thrust and Rudder Input 


Merion tag l = 


Card Columns 


1 is 


G25 


Zio 


4 an T= 60 


Option tag 2 = 


Format 


BLO 


EEC 


P10 


810 


SEO 


STB Thrust input 


bl aaa i 


(C0 iw Control tag and Option 
Vavue Or Constant thrust for the 
Starooard Screw if Variable thrust 


PorniOt Gesared 


Number of data points in the 


theust map. Left blank if 
constant thrust is used. 
Sreeethrust coerficient 

Data points for the independent 
variable, time in ascending order. 
Corresponding data points for 


the dependent variable. 


POttrememiaisi: Live uis 
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Cara Columns Format ey alte Tay 


I Las T 01602 - Control tag and option 
6-15 F10 As above for the port screw. 
16-25 F1O " i" " " " " 
26-35 F10 i " " " " " 
| a 1-80 8F10 " " " " i , 
io 0 8F10 ” " " " " " 


Option tag 3 = Rudder motion input 


card Columns Format Entry 
1 Joe i Hews e- COnLrol tag and option 
6-15 F10 VoOLVe Or snudaer lt constant rudder 


is desired 

Gs Fi0 Number of data points of the 
rudder map. Left blank if 
constant rudder is desired. 

Pe. i) SF10 Values of the independent variable 

time in ascending order. 

oi '5 0 8F10 Corresponding values of rudder angle 

| in degrees. Note: Positive rudder 


angle is right rudder. 


Brock 17 Bow and Stern Seal Pressure Differences Input 


Option tag 1 = Bow Seal 


Bard Columns Format Enver y 
al: Les i Ol7 0 Control tag and option 
6-15 F Number of data points for the 


bow seal 
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Gard Columns Format 
ee « 1-380 SEF10 
5 aa 1-80 8SF10 


Bie Y 
Time points in ascending order 
Corresponding differences in 


bow seal pressures 


Block 17 shall be left blank if a constant pressure is desired. 


miieae Constant shall be input in Block 5 and 6. 


However,- 


One seal may be held constant and the other varied. 


Option 2 = Stern seal 


Card Columns Format 
i 1-5 an 
6=15 F 
Pao. 1-80 SF1O 
Oy] s 1-80 SF10 
Block 18 Mithe Gard 


No option tags 


Card Columns Format 
iL 1-5 iL 
2 1-80 20A4 
mlock 19 %Fan Maps 


Option tag = 1 


Emery. 

Oly O0Z)=—Conerol tag and Option 
Number of data points for the 
stern seal 

Neue Ope Oils, Hon Che stem 
seal 

Micron Operon 1, f£or the stern 


seal 


Entry 
Gas O0)— Control tag 
Alphanumeric data to be printed 


as title in input summary 


Bow Seal Fans 
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Card Columns Format 
Ih 1-5 L 
6-15 I 
16-25 EF 
26-35 F 
36-45 F 
eee. 1-80 So 
After 1-80 8F10 
Past 
pressure 
data 
eard 


Mpe1on tag = 2 


Cousamiys 


Card Pommat 
iL 125. T 
-6-15 F 
LGD F 
26-35 EF 
ee. 1-80 8F10 
AEter 1-80 8F10 
last 
pressure 
Gata 
card 


Option tag = 3 


Entry 

Ole control tag and option 
Number of bow seal fans 

Speed of bow seal fans, rpm 

Number of data points (maximum of 25) 
If this is a 1.0 then the following 
Gata points are read. If left blank 
no further read takes place. 

Tabular value of fan pressure dif- 
ferential, psf. 


Gomereecoonai ug volumetric flow rate, 


GU vtt.7 Sec. 


Cushion Fans 


01902 - control tag and option 
Number of cushion fans 

Speed of cushion fans, rpm 

As option 1 above 

Tabular values of fan pressure 
differential, psf 


Corresponding volumetric flow rate. 


Stern Seal Fans 
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eard Cou amraays Format Entry 


i a) JL O190S  -Seemtrol tag and option 
6-25 F Number of stern seal fans 
16-25 F Speed of stern seal fans, rpm 
26-35 F Number of data points (maximum of 25) 
36-45 EF As option 1 above 
2 oa 310 8F10 Tabular values of fan pressure 


Gipaterential, ps t 


After [hess 8F10 Corresponding volumetric flow rates 
last 

pressure CUiasiin./ SCC. 

data 

@eard 


Note: The tabular pressure data for all the fans must be 


Stored in order of increasing pressure. (See discussion in 
manction FGl). The data is for a Single fan. 
Pemem@ecture Of Input Deck - The input deck for a given case or 


series of cases may be considered as being divided into two 
parts, "basic" @atemwana seeciite Gata. Basic data is defined 
as those inputs which are required to adequately represent 
the craft regardless of the type of case being considered, 
Data such as mass properties, craft geometry, fan maps, etc. 
fall into this category. Specific data is defined as those 
inputs which are needed to define a specific problem such as 
runs in waves vs. calm water, rudder motion vs. no rudder 
moOclLOn, etc. 

idee hecdtas which constitute the basic data are 


summarized below by control and option tag along with a 
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Meer description of the input. An asterisk indicates that 


more than one input card is read in under ENe wah chucomtre |. 


eae Option tag. 
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Required Input Summary (Listed in Numerical Order) 


Control Tag and Option WeSeCrlptELon 
OOds0: Timing information 
00103 Number of integrators and Ciro 


fOUerance s 


vOZ0) Vor Mass and inertias or mass dstribution 
00202* 

00300 Craft geometric data 

00401 Sidewall basic data 

00500 stern seal basic data 

00600 Bow seal basic data 

00701 Plenum geometry 

00702 Froude number at hump speed 
00801 Propulsion system basic data 
0090] | - Rudder geometry 

01990 . Aerodynamic reference lengths 
01200 iniettal conditions 

01300 End of data for one case 
01400 ' Fnd of all cases 

nO L600 (EME E ale 1 olebne 

01700 Rudder motions input 

01901* Bow seal fan maps 

01902* Cushion fan maps 


01903* Stern seal fan maps 
a ee 
*More than one input card under this control tag and 


ODEON . 
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Optional Input Summary (Listed in Numerical Order) 


Bontrol Tag 


Pee gig iim tiLOn lr Ayo t 


and Option Description Card Is Not Used 
pot02~* Subroutine NOUpreinerole FEOM Mndiyvrdual 
print switches subroutine 
00104* Tabular summary No tabular summary 
print switches 
00105 Option swithces Six Degrees of freedom, fixed 
for lateral plane thrust (uniess engine out is 
motions only, requested, surge velocity 
trim calculations allowed to vary 
and constant speed 
runs 
ot Ol* Wave input Calm water 
01800 ibavied Liowmrersl ats! No title printed 
feta pie Case Runs - In running multiple cases it 1s important 


that the user remember that any input data that is not read in 


again will be repeated from the previous case. Consider the 


Case where a run is to be made consisting of a run in waves 


followed by a turn in calm water at a different speed. Assume 


both intermediate and tabular ouput is desired for both cases. 


in terms of control tags the deck might look as follows: 


00101 
WOT 3 
WOZ0 2% 
00401 
00500 


00600 


Required Input 


*More than one input card under this control tag and option. 
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00701 


00801 

O20 1 

01000 

01200 

O16 00% 

01700 

SIIEMG Wes 

UENO es 

OARS MO) Shes 

ONGIINIO FA: Intermediate print 
00104* Lapultar print 

ClO 2= Wave input 

01300 End of case 1 
01200 ley sini teale Conditions 
01603 Rudecer MorlLon 

CaO Ze Wave input 

OA. 010 End of case 2 
01400 Bnew or fun 


mies purpose of the wave input card in case 2 is to read in 
the number of waves aS zero (calm water). Without this card 
case 2 would become a turn in waves. Similarly, if a run 
were to be made with case 1 being an engine out run and case 
2 being a turn with both engines operating, case 2 would 
have to reset the engine out time to a time larger than the 


run finish time by reading a large time on an 00802 card 





*More than one input card under this control tag. 
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(if case 2 had a finish time less than the engine out time 


from case 1 it would not be necessary to read this card). 


Output - The output from INCON is printed after an end of 
mse Card (01300) is read. It is a combination of input 
data and calculated results. A sample case is given in the 


Miput at the end of the report. 
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Subroutine INTGRL 


Discussion - Subroutine INTGRL is used to integrate a system 
Memerrest Order Ordinary differential equations. It uses a 

Variable time step technique based on the Runge-Kutta-Merson 
algorithm. A fixed time step R-K also available - See Block 
meoption tag 3 for 5 discussion on its use. The program will 


Stop the calculation if the time step becomes smaller than 


107° Sec. 
Usage |= CALL INTGRL (TIME) 


Descriptions of Parameters 


TIME = current value of the independent variable 


mep@ee- the input to subroutine INTGRL is controlled by the 


00103 control tag as discussed in block 1 of subroutine INCON. 


Output - If the integration step size is halved, INTGRL will 
mrant the current voluce On tame, the new Value of the time 
step, the integrator number, the calculated error and the 
input..maximum error as shown in the fame lewOueplr at the vend 
Smee report. This output is useful as it tells the user 
which integrator is causing the step size to be reduced and 
Meayeindicate an excessively tight input error tolerance 

fee discussion of 00103 card in discussion of subroutine 


INCON) . 
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Subroutine PROP 


Miecussion - Subroutine PROP is used to calculate forces and 
moments on the craft due to the propulsion system. The sub- 
moutine has the capability of uneven thrusting and receives 
jes information by using the Function Program FG1 to linearly 


Eiterpolate a map input in block 16. 
moage — CALL PROP 


Input - The input to subroutine PROP is controlled by block 


meee block 16 of subroutine INCON and is discussed there. 


Seeput —- The output from PROP is controlled by print option 
Meee LPROP. If this switch is 1, the program will print 
the calculated forces and moments as shown in the sample 


Output at the end of the report. 
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Subroutine RHS 


Mescussion - Subroutine RHS is the subprogram containing the 
FORTRAN expressions for the right hand side of the system 

of first order differential equations. It contains the logic 
to calculate the bubble volume, area and pressure as a function 
of time as well as the logic which sums the individual forces, 
moments and leakages from the various craft components to get 
Meet tOtals. RHS also contains the statements used for the 


feeeeing Of the scratch file for the tabular output. 
Usage - CALL RHS (VALUE) 


Mee iotion On aasaecers. 
VALUE = The array containing the values of the right 
Weancdeside Of the ditferential equations as 
Salemlatoawin  RHoG yy 1.62, hor a System Of 
Sie eee oor eenc efor yaar (y4,t)) the VALUE 


ede Weece Lomstora the values of f. (yj rt) 


faput — None 


~~ 


Suteut — The output from RHS is controlled by print option 
Byatehn THRS. When this switch is 1, the subroutine will 
print the output shown in the sample output in the back of 


mae report. 
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Simo routine RUDDER 


Discussion - Subroutine RUDDER contains the logic for cal- 
culating rudder forces and moments as well as rudder motions. 
Rudder motions are calculated by employing the Function Pro- 


gram FGl to linearly interpolate a map input in block 16 of 


INCON. 
Usage - CALL RUDDER 


Sout —- The input to RUDDER is controlled by block 9 and 16 


in Subroutine INCON and is discussed there. 


Poems — The output from rudder is controlled by print option 
Meee tRUD. If this switch is 1, the program will print out 
mies total forces and moments due to the rudder as shown in 


the sample output iaeene back of the report. 


Sh 





Subroutine SIDEWL 


Meeussion - Subroutine SIDEWL calcualtes the forces and 
moments acting on the craft due to the sidewalls as well 
as the leakage flow rates associated with any gaps which 
Open under the sidewalls due either to craft motion or 


Paves. 
wSeage - CALL SIDEWL 


mput —- The input to subroutine SIDEWL is controlled by 


blocks 3 and 4 of subroutine INCON and is discussed there. 


Peeput..- The output from SIDEWL is controlled by the print 
Option switch ISIDWL. When this switch is equal to 1, the 
program will print the sidewall gaps, immersion depths and 
total ar moments due to the sidewall-appendage 
combination as shown in the Cutpmi listwmg at the back of 


the report. 
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Subroutine STNSL 


Me cussion - Subroutine STNSL calculates the forces and 
moments acting on the craft due to the stern seal as well 
as the leakage flow rates arising from any gaps whicn open 
under the seal. The craft equilibrium leakage area is 


meemmed tO be located in the stern seal. 
Meage —- CALL STNSL 


S@eome —- The input data for STNSL is controlled by blocks 3, 


9, and 17 of subroutine INCON and is discussed there. 


Seeene — The output from STNSL is controlled by print option 
Switch ISTNSL. When this switch is equal to 1, the program 
Will print the seal wetted lengths, and gaps as a function of 
Station as well as the craft forces and moments due to the 
stern seal. An example is contained in the sample output at 


Mme back of the report. 


lps, 





menctions Tl and T2 


Discussion - Functions Tl and T2 are two small Function sub- 
programs used to calculate various trignometric relations used 


by subroutine WAVES. 
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Subroutine WAVES 


Mieeussion seoueroutine WAVES calculates the wave forces 
and moments acting on the craft. It also generates the 
wave heights at the various Stations around the seals and 
Sidewalls, as well as the bubble volume lost due to wave ele- 
vation. 

The sea state for irregular Seas is computed in WAVES 
by adding together a series of regular waves with an ap- 


propriate distribution of amplitude and frequency. 
Usage ~ CALL WAVES 


Input ~ The input to subroutine WAVES 1s controlled by block 


11 of subroutine INCON and 1S discussed there. 


Output - The output from WAVES is controlled by Pei ite Opiate 
Switch IWAVES. When this switch is set equal to 1, WAVES 
will print the wave elevation around Gee Sea! ceencdesigevalls, 
the wave é6levation at the center of Gravity and the volume 
meeuction of the pubis plenum due to waves as well as the 


total forces and moments acting on the craft due to waves. 
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USERS MANUAL FOR THE PLOTTING 
PACKAGE 


(INCLUDES THE MULTICURVE OPTION) 


maoCRIPTION: 


ies plotting package of the COLFIL subroutine is a 


meeehole routine that allows the user to: 


Ly 
2) 


3) 


7) 


Plot on either the CALCOMP or the printer 

Plot up to 3 curves per graph on the CALCOMP 
Intermix multicurve CALCOMP and single curve 
CALCOMP plots 

Revise the variables plotted on SuccesSive runs 
Choose from up to 26 different variables 

Print up to 16 variables per run 


Revise the variables printed on successive runs 


The method of employing this routine must be adhered 


emer 2CYrOrs will result. JCL for calling this routine 


Mm —rcnnd at the end of this manual. 


MEANING OF THE CARD 


1. Block 20 


Card Columns Format fe easy, 
1 1-3 iS 020 
2 iieemcacc muse be used 1f card 1 appears. see 


ample input 


1-10 Il HNeUPy  aeethe number Of this 


curve on the corresponding 
Capi len eco. | pertains to 


Capi mecca 2 — Graph 2, etc. 


EGiZ 





Card Columns Houma t me cy 

A “O“ must appear in these 

CO uncer Only One curve per 
GudiPieeaeca be plotted. "1" mact 
AOpen Gets MSsthestrirst oF 
MOGs SElvalimOnercurye per Graph, 
oe eee cets the anctermeditate 
aye wanideam o> if Ehls 1S the 


Pao ecwrve  onecias Graph. 


fo Block 21 


Card Columns Format Pinay, 
1 8: Ls O21 - Causes the summary variables 
to be read 
2 1-16 812 Punch in the numbers corresponding 
to the variable names desired to 
be printed in summary one. 
3 1-16 Sal The same aS card one but, for 
Summary Two, these cards need 
not be repeated for successive 
runs unless changes are desired. 
ee Block 22 
Card Columns Format Petters Ys 
i 1-3 ES 022 - Read in the number of grapns 
desired, whether or not Calcomp is 
Gdectieeceeeread the X and Y axes for 


the graphs and read the first line 


EGS 





card 


e@ kumins 


Format 


li 


2 


22 


6A8 


Entry 
of the graph title. 

The number of graphs desired right 
justentCd, lp tO and including 10: 
A zero if print plot is desired or 
ewer Shite justi tLed, 1f£ CALCOMP 

is desired. 

The numerical values of the X and 

Y axes alternating in the following 
Manner. The X axis of the first 
graph then the Y axis of the first 
Fae pimaiG —COMLInuIng for all 
desired graphs. A table relating 
axis number and name is at the 
endorses, paler. 

The Alpha-numeric title that will 


appear on the first line of the 


SG Leo slots [it CALCOMP is not 


desired, replace this card with a 
blank card. Your name should appear 
somewhere on this line for computer 


identification purposes. 
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020 


O22 
0201 .——____________ 
ewes 
eeeeoucput for LEO 
Standard input deck 
013 


020 
Z2 


O21 Gye 


BOLG 
piel5 


Changes desired 


Changes desired 


013 
O14 


Se eee lei Oe DECK 


tiewtiteibeates that this 1S the first of 


multiple curves on graphs l and 2. 


swat yeOnc wihivorant the Linst two vartables, 
Oi, 15, and summary two the following two, 
Cie Gr 


Ivo Gilapmicwapescalled £0r, and a CALCOMP output 
CusphaeieploessO!, 15, Time and X displacement 


Grape 2 plocs 01, 16, Time and Y displacement 


As Above. These can be changed on subsequent 


runs. 


As Above. 
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RELATION OF VARIABLE NAME AND NUMBER 


NUMBER 


Ol 
02 
03 
04 
05 
06 
07 
08 
09 
10 
ea 
a2 
IL 
14 
15 
16 
Ly 
(18 
LS 
20 
Zi: 


22 


NAME 





Time 

Wave Height 

Z Displacement 

Patch, Angle 

Plenum Pressure 

Bow Acceleration 
GrG. -Acce leratyon 
Fan Power 

Roll Angle 

Yaw Angle 

Lateral Acceleration 
Speed through the water 
Turn Radius 

dMGWoigl) Lesher! 

X Displacement 

Y Displacement 

Air Flow in 

Pane 2 low vOut 

Net Force X Direction 
Wave Force X Direction 
Net Force Y Direction 


Net Force Z Direction 
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NUMBER NAME 





23 Net Torque X Axis 
24 Net Torque Y Axis 
25 Net Torque Z Axis 
26 Rudder Angle 


We 





APPENDIX A-A 


CHANGING THE MODEL 
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Ehengang the Model 
Appendix A-A 

meroduction — 

AS previously stated, the program has been written 
in a modular form to facilitate changes to the program 
which reflect major design changes to the craft, i.e., 
a Gifferent type of seal, addition of controls, etc. Since 
it is obviously impossible to program all possible con- 
figurations in one program, these types of changes will 
require either replacing a whole module (sSub-program) or 
altering the Fortran statements within a given sub-program. 
In order to achieve this, a more detailed knowledge of the 
structure of the program and the method of data trans- 
ference between subprograms is required. This appendix 
will attempt to give a brief explanation of the program 
Structure and will then discuss the details of how one 
wollild make certain anticipated changes. 
General Structure 

The primary variables (displacements, velocitites, pres- 
Sure, etc.) are contained in an array of length 40 words, 
named VAL; which is transmitted through the various sub- 
programs using a labelled common block called VARBLE. Within 
each subprogram the various elements of VAL are equivalenced 
to local variables, which are given mnemonic names which 
closely resemble the standard mathematical symbols used for 


the variable (e.g. u, v, w, Ds become U,V, W, PB in the 


IOS, 





Fortran listing). The second through the sixteenth elements 
Seve, are reserved for the output of integrators and must 
be used in numerical order, i.e., when adding a new equation 
fo be integrated, the output of the integrator must be equi- 
femenecead tO the first unused location in VAL (VAL(12) for 
mie Current program). 

wie force and moment components acting on the craft are 
calculated for each element in the appropriate subprogram 
and are transmitted through various labelled common blocks. 
These force and moment components are then summed in sub- 
routine RHS to get the appropriate acceleration and velocity 
components. The accelerations and velocities are trans- 
mitted back to subroutine INTGRL in an array called VALUE. 
HNTGRL then integrates these first derivatives to get the 
new velocities and displacements. 
Anticipated Changes 
eee placing Analytic Coefficients or Curve Fits With 

Tabular Experimental Data 

As experimental data is obtained it may be desired to 
Peelace coefficients based on analytic calculations or 
Seeiations from Simple shapes by the experimental data. 
This data will often be tabular in nature. For functions 
of a Single variable this is quite simple to do as the 
necessary table look-up routine already exists within the 
B2ogram (function FG1). The user would be required to 


make the following changes. 
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[ee oeaeements LO read in the tabular data must 
be added to subroutine INCON. 

Zee common block must be provided to transfer 
the data to the subroutine where it is to be 
used. 

3. A statement must be added to call FG1 to do 
Pie weal sLOOK—Up. 

A good example of the above procedure is afforded the 
Meer DY the fan maps as currently programmed. This logic 
uses tabular data and the user should refer to input block 
19 of subroutine INCON and subroutine FAN to see an example 
of the three steps outlined above. 

b. Adding a Stern Seal Equation of Motion 

Pee wcurrent program has a stern seal model which is 
based on the assumption that the seal is completely de- 
formable ine sthimenese. and as such there is no dynamic 
equation of motion. For seal designs where there is an ap- 
Mraiciable degree of stiffness a seal equation of motion 
should be added. 

Assume that the seal equation of motion to be added is 
a standard spring-mass-damper syStem and therefore may be 
represented by the following system of first order dif- 


ferential equations. 


Gr. = (-d@, - kO, + M.)/m (1) 
ee Ob. (2) 


aby lk 


The following program changes would have to be made in 


Seaer to incorporate Equations (1) - (2) into the current 


program. 


de 


Pegie must be addea to block 5 of subroutine 
INCON to read in and store the constants da, 
Mmancieas wellewas the initial conditions 
Lor 9. and Q.- Care should be taken in the 
case of the latter to eave these initial 
conditions using Fortran variables with names 
other than the names which will contain the 
current (running) values of the seal displa- 
cement and velocity. This is due to the fact 
that when running multiple cases the program 
will use the currently stored values of 8. and 


0. PEeSimmeneeGladmecase as Initial conditions for 


BiemNoweeaceclnbessetnie initial conditions are 


reset. An example of this can be seen on page 
ioc maEogrameltSting akter statement num- 
Pomme WOswneremthe Crake initial conditions 
Meemacad mime Saving the anitial conditions 
under a different name also has the advantage 
of allowing the input to be in degrees rather 
than in radians aS required by the program, as 
the conversion can be done ine arithmetic 


replacement statements added in step 2 below. 


eZ 


iMimebeck 1S %ef INCON, (after statement 1300 
but before statement 1400) arithmetic re- 
placement statements should be added to con- 
eG Chic tlal condi tlonsato radians and 
radians/sec. and then store these values in 
the running variables for 8, and 6. (see page 
Z05 of thie program listing after statement 1302 
for an example). Note that in the program 

the variable RAD is the number of degrees in 

a radian and may be used for the conversion 

of units (see example cited above.) 

The variable names for the seal velocity and 
displacement should be equivalenced to VAL (12) 
am@e VAI (13) respectively. This automatically 


makes these variables the output of integrators. 


“This equivalence statement should appear in sub- 


routines INCON, RHS and STNSL. 

The number of differential equations on the in- 
put card after the 00103 card should be increased 
Cone mEeel?, sandicating that two additional 
Qifferential equations have been added. In 
ad@atiom, error tolerances for these integrators 
must be given (see pagel54 of this report for 
emomscucsstonvor this point.) 

_eGormmenmo lockmror transmit d, Kk, and m must be 


added to subroutines INCON and RHS. 


is 





Geel aesubroutine RHS the Fortran equivalent of 
the right hand side of Equation (1) must be 
set equal to VALUE (11) while on must be set 
eollakeuo VALUE (12). (See examples in listing 
Of Subroutine RHS.) 

7. In subroutine STNSL the expressions for Mos the 
total applied moment on the stern seal aout its 
hinge point must be added. These are placed in 
STNSL as this subroutine has the logic for cal- 
culating wave elevations relative to the seal, 

SEC lin acoltlon,. Lhe equations for the stern seal 
leakage gap is STNSL must be modified to take into 
account the seal motion. 

8. A common block must be added to subroutines STNSL 
an RHS to transmit the value of the moment Mo to 
RHS. 

c. Adding New Forces oe MoMCMeseACElNG On a Crart 

If the modification of the program entails the addition 
Sreany new force or moment component acting on the craft, care 
must - Eig te eomacc sents COMmpOnentato the bending moment cal- 
culations as well as the craft motion calculations. If this is 
not done the equilibrium check done by subroutine SAM prior to 
calculating the bending moments slate balance and the cal- 
culation will be terminated. 


The following changes would have to be made: 
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In subroutine RHS the new force component should 
be added to the WRITE statement preceding state- 
ment number 111 on page 218 of the program listing 
as well as into the calculation of GF in RHS. If 
the force component has a fixed moment arm, a 
common block should be added to transmit its value 
to subroutine SAM. If the moment arm is time va- 
rying the moment arm should also be added to the 
WRITE statement. 

In subroutine SAM the force and arm (if necessary) 
should be added to the following READ or WRITE 
statements; page 220, statement numbers 42 and 


44; page 220, after statement number 6. 


The force and moment components must be added 


EGethe Cquilibrium check on the top of page 


The force and moment components must be added 

Pomene calculation o£ etme Jongitudinal and 
lateral moment and shear calculations. In 

these calculations the variables XMI(IX) and 
VWiktG@hy vane ethe drStance to the current plane 

See erest see tie Uses should Study the current 
mEeeciane listing (pages L863 —- 233 to ascertain, ene 
procedure used to add these force and moment com- 


WOneiito 


GS 


~ ae 





Jee) estoy NB ILD,¢ /8\ealB: 


ABNORMAL TERMINATION MESSAGES 
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APE END EX A-C 


SIDEWALL INTEGRAL TABLE PROGRAM 


DiSCuSo ION 
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Sidewall Integral Table Program 


Appendix A-C 


The sidewall Integral Table Program evaluates tables of 
integrals of sectional hydrostatic and hydrodynamic pro- 
perties over the whole and partial iengths of a SES side- 
fer. These tables are calculated as a function of sidewall 
Meet, trim and length of a regular wave. The tables are 
then subsequently accessed by the main SES Motions and Load 
Program and when combined with knowledge of both the craft 
and waves instantaneous position and motion produces the 
hydrostatic and hydrodynamic forces and moments acting on 
the sidewalls. 

The primary input to this program is a table of water- 
lines and pees roneinc beams for each chosen station with 
arbitrary spacing in both waterlines and stations. The pro- 
gram is thus applicable to a general sidewall form. In 
addition to this input, the user must specify the range of 
drafts, trim (pitch) angles, and wavelengths (encompassing 
all Pereonents of a specified sea state) that the craft is 
expected to range thru in the subsequent simulations per- 
formed by SES Motions and Loads Program. The wave-lengths 
are specified via the wavelengths parameter (27I/A). The 
first value of this parameter must be zero which corresponds 
to an infinite wavelength and is used to SIDEWL subroutine 
of the SES Motions and Loads Program to determine hydro- 


Static and hydrodynamic forces due to body motion alone. 


ds ss 





the tables at the remaining values are used by the WAVES 
Subroutine for wave induced forces. Since the center of 
Gravity is not necesSarily known at the time thie program 
1S executed, the user must also specify a BONG ieee aves! 
reference coordinate measured from the transom about which 
Ene integrals are to be evaluated (a value near mid-ships 
1S recommended). The range of drafts and trims at the longi- 
tudinal reference point must be selected in such a manner 
So that the top of the plenum does not sink below the water 
Surface. 

The description and order of Ehemcards of tEhe input 
data deck is as follows: 


Best Card 


Columns Format Entry 
5 . iz | No. Of wave patameters 420 
ot a F increment in wave parameters (ft.71) 
LSE) EF | first value of wave parameter (ft. 7+) 
16-20 i new Of drafts <5 
21-25 : | increment in draft ‘Cael ae 
ie 30 F first value of draft (in.) 
om | il no. of trims <7 
56-40 F InecEeKeCmenin trim (deq. ) 
41-45 F Eaves vane Ot trim (deg. ) 
46-50 ai ho. of stern beams 436 
51-55 F increment of stern beam (ft.) 
56-60 F first value of stern beam (ft.) 


180 





ema Card 


Columns 


lame, 
SLO 


Sea Cara 


Solumns 


1-5 
Seah Card 
Columns 


1-10 


Siem Card 


Columns 








15 


Next NSTA Cards (one card 


Columns 


l= 5) 


Format 


EOrpmar 


ig 


Format 


Format 


Ors 


Oma 
F5 
eS 
F5 
FS 


F5 


HEME ey 
NOeeOreStatlons = NSTA £95 


Wee OrmOl FGabhis pcr Station <8 


Peasy, 


distance from keel to top of plenum (ft.) 


ENG ry 
distance of reference position 


freomveransom. (it ) 


BNEry 

New -OtimimEecdral. limits sin addition 
Bemete it SCN Ovets he whole eng tly ai 
integral limits (distance from transom) 
(etter) 

per station) 

Eiieay 

distance of station from bow (in.) 
lowest waterline beam (ft.) 
corresponding height of beam (ft.) 

next waterline beam (ft.) 


corresponding height of beam (ft.) 
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Columns Format Entry 


66-70 F5 iNSSeEe Wiehsicvellsiavel Jerevcii (Geeme) 


7 lls F5 corresponding height of beam (ft.) 


This program transmits the calculated tables to the SES 
Motions and Loads Program by writing the tables on units which 
Meee later read by the main program. Unit number 10 contains 
the tables obtained by integrating over the entire length and 
Mme succeeding units contain the tables obtained by integrating 


Over partial lengths. 
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APPENDIX B 


FORTRAN PROGRAM LISTING 
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